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leakage,  face  temperature,  and  wear  are  measured. 
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APPLICATION  OF  RESEARCH  TO  THE  NEEDS 
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Mechanical  face  seals  are  used  in  numerous  applications  in 
Naval  machinery.  These  applications  range  from  propeller  shaft 
seals  to  boiler  feed  pump  seals.  In  such  equipment  the  mechanical 
seal  plays  a vital  role.  When  such  seals  fail,  repair  is  costly 
both  in  terms  of  lost  time  and  direct  cost,  so  any  improvement  in 
seal  life  and  reliability  would  be  of  significant  benefit. 

As  more  advanced  equipment  is  designed,  it  is  sometimes  diffi- 
cult to  achieve  desired  performance  in  more  severe  service  environ- 
ments with  the  present  state  of  the  art  of  seal  design.  Thus,  an 
improvement  in  seal  technology  would  serve  this  important  application. 

The  immediate  objective  of  the  research  herein  is  to  further 
the  understanding  of  mechanical  face  seal  lubrication  phenomena. 

The  ultimate  objective  is  to  develop  the  capability  of  designing 
contacting  face  seals  having  a longer  life,  greater  reliability, 
and  for  extreme  environments.  Thus,  the  objectives  of  this  research 
are  compatible  with  mechanical  face  seal  needs  for  Naval  machinery. 


CHAPTER  1 


LUBRICATION  IN  MECHANICAL  FACE  SEALS 


INTRODUCTION 

In  applications  where  a rotating  shaft  must  pass  from  one  fluid 
region  to  another,  contacting  mechanical  face  seals*  play  the  essen- 
tial role  of  minimizing  the  transfer  of  fluid  between  the  regions. 
Applications  of  face  seals  range  from  water  pump  seals  to  process 
pump  seals  to  propeller  shaft  seals. 

The  performance  and  reliability  of  contacting  mechanical  face 
seals  are  of  great  importance  for  any  type  of  equipment  where  minimal 
leakage,  high  reliability,  and  long  life  are  necessary.  Even  for 
equipment  where  these  factors  are  not  so  critical,  seal  failures  and 
short  seal  life  lead  to  high  operating  cost  due  to  down  time  and 
maintenance  cost. 

Even  though  mechanical  face  seal  technology  has  been  steadily  im- 
proving over  the  past  several  decades,  further  improvement  in  the  state 
of  the  art  of  seal  design  would  be  most  beneficial.  Although  seals 
having  an  acceptable  life  and  reliability  can  be  designed  for  many  ap- 
plications, further  improvement  in  seal  life  and  reliability  would  re- 
sult in  significant  cost  savings  to  the  user.  Also,  there  are  numerous 
mechanical  face  seal  applications  where  seal  loading,  reliability,  life, 
and  leakage  requirements  are  difficult  to  achieve  within  the  present 
state  of  the  art.  Examples  of  such  applications  are  seals  for  pumps 

*The  class  of  low  leakage  face  seals  where  there  is  definite  contact 
and  wear  of  the  faces  as  opposed  to  hydrostatic  or  hydrodynamic  where 
a definite  clearance  is  maintained. 


for  nuclear  power  plants  and  seals  for  large  diameter  submarine  pro- 
peller shafts.  Additionally,  the  friction  losses  in  face  seals  repre- 
sent a significant  fraction  of  energy  consumed  for  pumping  purposes, 
within  the  present  state  of  the  art  of  seal  design,  it  is  very  diffi- 
cult to  design  a low  leakage  seal  that  also  has  a low  friction  loss. 

The  main  barrier  to  the  advancement  of  the  state  of  the  art  is 
that  the  mechanics  of  seal  operation  are  not  well  enough  understood 
to  be  able  to  reasonably  anticipate  seal  performance  as  a function  of 
design  parameters.  There  are  no  well  established  fundamental  theo- 
retical bases  that  can  be  used  to  indicate  the  type  of  seal  design 
that  will  give  improved  performance . Improvements  that  have  been 
made  have  been  brought  about  largely  by  trial  and  error  combined  with 
elementary  sealing  theories. 

In  order  to  be  able  to  predict  the  performance  of  contacting 
face  seals  as  a function  of  design  parameters,  it  is  essential  that 
the  lubrication  mechanisms  between  the  faces  be  well  understood.  At 
present  it  is  known  that  hydrodynamic  or  hydrostatic  lubrication  plays 
some  role  in  providing  load  support  for  oil  seals  as  well  as  water 
seals.  But,  the  precise  nature  of  this  lubrication  is  not  known. 
Several  theories  have  been  put  forth.  However,  these  theories  have 
not  been  verified  for  contacting  face  seals,  and  it  is  not  possible 
as  yet  to  use  these  theories  for  the  design  of  contacting  face  seals. 

In  this  work,  the  results  from  further  research  into  hydrodynamic 
lubrication  mechanisms  in  face  3eals  is  reported.  Much  has  been 
learned  about  these  mechanisms,  and  there  is  some  promise  that  effec- 
tive use  of  hydrodynamic  mechanisms  can  be  made  to  improve  face  seal 
performance . 
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MECHANICAL  FACE  SEAL  LUBRICATION 


The  mechanical  face  seal  consists  basically  of  two  annular  rings 
which  rotate  relative  to  each  other  and  which  are  pressed  together  by 
spring  and  fluid  pressures  (see  Figure  1-1)  . The  surfaces  that  rub  to- 
gether are  generally  manufactured  as  flat  as  possible  initially  so  as 
to  minimize  leakage.  The  effective  gap  between  the  faces  is  ideally 
quite  small  (order  of  1 urn)  so  that  leakage  flow  across  the  faces  will 
be  quite  small.  The  difficulty  in  designing  a mechanical  seal  is  in 
maintaining  the  gap  at  a very  low  value  while  at  the  same  time  pro- 
viding a definite  lubricant  film  between  the  faces. 

The  load  that  must  be  supported  at  the  faces  of  a mechanical  seal 
is  due  primarily  to  loading  caused  by  the  sealed  pressure.  The  unit 
face  load  can  be  expressed  as  some  fraction  B of  the  sealed  pressure 
where  B can  be  made  greater  or  less  than  unity  by  geometry  selection. 
The  load  support  at  the  faces  is  derived  from  fluid  pressure  and  me- 
chanical pressure . If  the  fluid  pressure  at  the  faces  is  large  enough 
to  support  all  of  the  load,  then  there  will  be  no  contact  and  no  ad- 
hesive wear.*  If  none  of  the  load  is  supported  by  fluid  pressure, 
the  load  must  be  carried  by  mechanical  contact  and  the  wear  rate  will 
be  large . 

In  practice,  seals  often  operate  at  one  of  two  extremes.  At  one 
extreme,  a large  gap  will  be  created  by  hydrostatic  or  hydrodynamic 
effects  and  the  seal  will  leak  a large  amount.  At  the  opposite  ex- 
treme, the  gap  will  close  completely  and  leakage  will  be  very  low. 


* 

There  may  still  be  abrasive  wear  even  if  the  surfaces  do  not  touch. 


However,  the  beneficial  lubricant  film  will  be  lost, and  wear  and  heat 


generation  will  increase. 

Based  on  the  above,  it  can  be  concluded  that  an  effective  seal 
should  operate  between  these  two  extremes — having  both  adequate  lubri- 
cation and  low  leakage.  To  do  this  requires  that  any  fluid  pressure 
generation  mechanism  used  to  provide  lift  and  lubrication  to  the  seal 
must  be  very  carefully  controlled.  At  present,  in  commercial  contact- 
ing face  seals  this  lubrication  is  left  primarily  to  chance,  and  often 
such  seals  operate  at  one  of  the  extremes  mentioned.  Quite  commonly 
such  seals  will  operate  in  the  low  leakage  condition  where  wear  pro- 
gresses at  a definite  rate.  Such  seals  are  quite  satisfactory  for 
many  applications.  However,  when  pressure  or  speed  requirements  are 
increased,  then  this  mode  of  operation  may  lead  to  rapid  failure . 

In  order  to  be  able  to  design  into  the  seal  the  proper  balance 
between  lubrication  and  leakage,  the  origin  of  lubrication  mechanisms 
must  first  be  fully  understood.  A considerable  amount  of  research  has 
been  accomplished  in  this  area  over  the  years.  Various  theories  have 
emerged  and  have  been  described  in  the  previous  annual  report  [1]*. 
These  theories  fall  roughly  into  two  categories,  hydrostatic  and 
hydrodynamic.  In  this  research  program,  attention  has  been  directed 
toward  studying  the  effects  of  waviness  as  a source  of  hydrodynamic 
pressure  in  face  seals.  Particular  attention  has  been  focused  on 
the  effects  of  roughness  and  low  viscosity  applications  where  some 
touching  is  expected.  Waviness  was  selected  among  the  various 

* 

Numbers  in  brackets  refer  to  List  of  References  at  the  end  of  report. 


hydrodynamic  theories  because  it  is  controllable  in  seals  and  a 
better  understanding  could  lead  to  in®>roved  seal  designs. 

Sources  of  waviness  in  seals  have  been  described  in  detail 
previously  [1].  Waviness  may  be  produced  accidentally  as  an 
uncontrolled  variable  or  intentionally  in  a precisely  controlled 
manner.  A wavy  seal  ring  is  shown  in  Figure  1-2. 

WAVY  FACE  SEAL 

To  describe  the  wavy  seal  lubricating  mechanism  in  more  detail, 
consider  the  ring  shown  in  Figure  1-2  as  being  pressed  against  a 
flat  rotating  ring.  Then  the  problem  shown  in  Figure  1-3  results. 

The  seal  is  subjected  to  outside  and  inside  pressure.  The  wavy  face 
(shown)  is  fixed,  while  the  flat  face  rotates.  Given  the  two  waves, 
for  example,  there  will  be  two  points  of  minimum  film  thickness  be- 
tween the  faces  as  shown.  The  seal  may  actually  contact  in  regions 
near  these  points  as  shown.  Angular  rotation  of  the  mating  face 
causes  hydrodynamic  pressure  to  build  up  in  the  converging  regions 
as  the  point  of  minimum  film  thickness  is  approached.  In  the  diverg- 
ing regions,  cavitation  will  occur  as  the  pressure  attempts  to  be- 
come negative.  This  cavitation  is  essential  for  a net  lifting  force, 
and  for  an  outside  pressurized  seal,  the  cavity  will  appear  as  shown 
in  Figure  1-3.  The  fluid  will  flow  in  streamers  across  the  cavity, 
and  a full  film  will  again  develop  at  the  downstream  cavity  boundary. 

The  problem  to  be  solved  is  to  determine  the  fluid  and  mechan- 
ical pressure  distributions  for  the  given  configuration.  Given  the 
pressure  distribution,  then  load  support  and  leakage  can  be  calculated 
and  the  effect  of  various  parameters  can  be  studied. 
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ROTATION  OF  MATING  FACE 
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Figure  1-3.  Wavy  Seal  Geometry 
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The  hydrodynamic  pressure  distribution  for  this  problem  has 
been  solved  by  various  methods  by  Findlay  [2],  Pape  [3],  and 
Stanghan-Batch  and  Iny  [4]  for  perfectly  smooth  faces.  These  re- 
sults show  that  even  a small  waviness  produces  sufficient  load  sup- 
port for  complete  liftoff.  However,  heavily  loaded  or  low  viscosity 
seals  do  not  operate  with  full  hydrodynamic  load  support  and  com- 
plete separation.  A definite  wear  results.  Surfaces  of  such  seals 
contact  during  operation,  and  hydrodynamic  pressure  is  affected  by 
interactions  with  surface  roughness  (See  Figure  1-4) . Seals  of 
this  type  operate  in  a mixed  friction  regime.  Hydrodynamic  load 
support  due  to  waviness  continues  to  provide  a significant  fraction 
of  the  load  support,  but  asperity  contact  must  provide  the  balance. 
Thus  these  important  effects  must  be  included  in  a wavy  seal  lubri- 
cation model. 

In  the  first  annual  report  for  this  project,  reference  (5] , 
the  above  problem  was  solved  using  a one  dimensional  theory.  In 
the  second  annual  report,  [1],  the  much  more  complex  two  dimensional 
solution  to  the  above  problem  was  solved.  The  effects  of  waviness, 
roughness,  asperity  contact,  wear,  cavitation,  and  elastic  deflec- 
tion were  included  in  the  model.  Using  this  model  predictions  were 
made  for  the  relative  wear  rate,  friction,  and  leakage  as  a func- 


tion of  roughness,  waviness,  speed,  size,  pressure,  viscosity,  and 


material. 


A number  of  conclusions  were  reached  based  on  the  studies 
reported  in  reference  [1] . Those  having  impact  on  this  present 


work  are: 
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Figure  1-4.  Seal  Roughness 


1)  The  effects  of  roughness  on  hydrodynamic  lubrication  are 


not  completely  understood.  Certain  fundamental  questions 
remain  concerning  the  roughness  model  used. 

2)  As  to  the  potential  of  utilizing  hydrodynamic  effects  to 
advantage  by  design,  the  results  show  that  wear  rate  and 
friction  can  be  greatly  reduced  while  maintaining  leakage 
at  acceptable  levels. 

3)  While  a comparison  of  predicted  results  to  experimental 
results  given  in  the  literature  is  generally  good,  data 
contained  in  the  literature  are  incomplete,  so  more  com- 
plete experimented,  data  are  needed  for  comparison. 

The  present  year's  activities  have  been  directed  toward  pro- 
viding additional  information  on  these  important  topics.  In  Chap- 
ter 2 the  effect  of  roughness  on  lubrication  is  reviewed  and  several 
new  developments  are  presented.  This  particular  topic  is  being  in- 
vestigated by  a number  of  people,  and  certain  questions  have  yet  to 
be  answered.  As  a consequence,  a seal  model  based  upon  a particular 
roughness  model  can  become  suddenly  obsolete. 

In  Chapter  2,  an  attempt  is  made  to  resolve  certain  fundamental 
questions.  While  contributions  to  the  knowledge  were  made,  other 
questions  were  raised  in  doing  so.  The  question  of  roughness  effects 
remains  unsettled  and  a final  seal  model  must  await  these  conclusions 
In  spite  of  this  limitation,  a different  surface  roughness  inter- 
pretation has  been  made,  and  results  compare  favorably  in  some  re- 
spects to  those  from  the  previous  model  [1] . The  interpretation 
made  also  points  out  the  strong  need  for  additional  experimental 
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work  which  is  being  pursued  as  discussed  below. 

In  Chapter  3,  attention  is  directed  toward  developing  a 
rational  design  procedure  for  a wavy  mechanical  face  seal.  The 
parameter  studies  of  the  previous  report  [1]  have  been  used  as  a 
basis  for  this  work.  Additional  studies  have  been  made  and  a de- 
sign procedure  has  been  established.  Results  are  presented  which 
show  the  advantages  of  a wavy  seal  as  well  as  some  of  the  pitfalls 
of  implementing  such  a seal  design. 

In  the  previous  annual  report  [1] , the  design  of  the  main 
portion  of  an  experimental  test  apparatus  was  described.  During 
this  past  year  the  design  of  the  auxiliary  equipment  for  the  ma- 
chine was  completed  and  the  entire  machine  has  been  built  and  tested 
out.  Mechanical  seal  experiments  are  now  being  carried  out  using 
this  apparatus.  In  Chapter  4 a complete  description  and  documenta- 
tion of  the  test  apparatus  is  provided. 


CHAPTER  2 


REEVALUATION  OF  SURFACE  ROUGHNESS  EFFECTS 
IN  HYDRODYNAMIC  LUBRICATIONS 

In  the  previous  annual  report  [1] , discussion  was  presented  in- 
dicating that  the  precise  effects  of  roughness  on  lubrication  were 
not  clearly  understood  and  that  the  correct  method  of  mathematical 
treatment  of  the  problem  was  still  developing  and  a subject  of  con- 
tinuing debate.  Given  the  importance  of  roughness  effects  on  lubri- 
cation in  a wavy  seal  con  figuration  as  pointed  out  in  the  previous 
annual  report  [1],  it  was  decided  that  some  attempt  should  be  made 
to  resolve  the  question  as  to  the  validity  of  various  approaches, 
so  that  the  best  available  models  of  surface  roughness  effects 
could  be  incorporated  into  the  wavy  seal  model. 

After  reviewing  the  latest  literature , it  was  decided  that  the 
recent  work  of  Patir  and  Cheng  [6]  held  the  greatest  promise  of 
providing  a useful  yet  accurate  mathematical  description  of  surface 
roughness.  Therefore  in  this  chapter  a comprehensive  review  of 
surface  roughness  theory  is  presented.  The  Patir  Cheng  [6]  work  is 
then  evaluated  using  the  theoretical  basis  established.  Using  tech- 
niques developed  herein,  the  limitations  and  uncertainties  concerning 
the  Patir  Cheng  work  are  developed.  Finally  some  conclusions  are 
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made  as  to  the  significance  of  these  results  on  the  rough  wavy  seal 
problem. 

As  a first  step  to  place  the  recent  work  in  perspective,  a 
literature  survey  is  presented. 


LITERATURE  SURVEY 

The  effect  of  surface  roughness  on  the  lubrication  of  hydro- 
dynamic  bearings  has  been  the  subject  of  research  only  in  the  past 
decade.  Tzeng  and  Saibel  [7]  developed  the  concepts  of  stochastic 
modeling  as  applied  to  a slider  bearing  possessing  transverse  rough- 
ness in  1967.  During  the  next  several  years  Christensen  and  Tonder 
derived  average  Reynolds  equations  governing  the  expected  pressure 
in  bearings  having  transverse  and  longitudinal  roughness  [8],  [9]. 

In  all  of  this  work  classical  lubrication  theory  was  applied  to 
rough  surface  situations.  Based  on  the  experimental  work  of  Ostvik 
and  Christensen  [10]  it  was  observed  that  after  some  run-in  and 
wear  in  roller  bearings,  the  surfaces  often  exhibited  a two-dimen- 
sional longitudinal  roughness.  Hence  the  Chris tensen-Tonder  work 
was  oriented  toward  surfaces  with  this  type  of  roughness  having  an 
approximate  Gaussian  frequency  distribution.  Ultimately  the 
Christensen  approach  was  extended  to  isotropic  surfaces  [11]  and 
this  treatment  is  considered  mathematically  invalid  by  Elrod  [12] . 

In  another  publication,  Elrod  [13]  points  out  that  the 
classical  Reynolds  equation  may  not  be  applicable  in  the  analysis 
of  surface  roughness  effects  due  to  dimensional  considerations. 

Elrod  coined  the  terms  "Reynolds  roughness"  and  "Stokes  roughness" 
to  emphasize  the  restrictions  on  the  classical  Reynolds  equation. 
Second  order  effects  of  the  Stokes  equations  to  account  for  large 
surface  curvature  were  included  in  an  analysis  by  Sun  and  Chen 
[14],  but  the  method  is  mathematically  laborious.  An  alternative 
to  the  use  of  Stokes  equations  is  the  development  of  a modified 
Reynolds  equation  which  includes  rough  surface  functions  explicitly. 
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Both  Elrod  [12]  and  Godet  [15]  have  derived  such  equations.  It  is 


emphasized  by  Godet  that  the  total  surface  roughness  may  not 
necessarily  be  lumped  onto  one  surface  as  was  done  by  Christensen 
in  accord  with  the  classical  Reynolds  equation. 

The  special  Reynolds  equation  was  used  to  derive  an  average 
equation  in  the  spirit  of  the  Christensen  equations  by  Cheng  and 
Patir  [6 1 and  by  Elrod  [16]  . Elrod  used  a perturbation  technique 
applied  to  the  modified  Reynolds  equation  which  solves  for  both 
the  average  pressure  and  the  standard  deviation  of  the  pressure 
value.  However,  the  result  is  not  applicable  to  bearings  where 
contact  occurs,  and  an  estimate  of  the  error  for  closely  approach- 
ing surfaces  is  not  made.  The  average  flow  model  developed  by 
Cheng  and  Patir  allows  for  contact  of  the  bearing  surfaces  and 
is  applicable  to  real  three-dimensional  surfaces.  It  was  found 
that  the  average  flow  model  yields  consistent  results  in  the  limit 
with  the  Christensen  two-dimensional  roughness  work.  Gaussian 
surfaces  were  used  for  these  results  and  the  autocorrelation  and 
non-isotropic  character  of  the  surfaces  were  imposed  by  use  of  a 
method  based  on  the  joint  ideas  of  Peklenik  [17]  and  Patir  [18]. 

The  work  of  Cheng  and  Patir  also  discredited  the 
Christensen-Tonder  isotropic  model  by  showing  that  the  latter 
results  underestimated  hydrodynamic  load  support.  Although  the 
longitudinal  model  of  Christensen  is  considered  a close  approxi- 
mation to  real  surfaces , research  by  Lebeck  [ 1 ] on  mechanical 
face  seals  points  out  that  since  side  leakage  is  prevented  in 
contact  regions,  the  Christensen  longitudinal  model  predicts  un- 
realistically large  values  of  pressure  when  contact  occurs. 
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While  the  preferred  treatment  of  surface  roughness  in  hydro- 
dynamic  systems  when  contact  may  occur  is  definitely  not  decided, 
the  i-atir  Cheng  Method  does  hold  considerable  promise.  However, 
results  from  the  average  flow  model  are  presented  as  a unique  solu- 
tion and  the  user  is  not  given  warning  about  the  limitations,  ac- 
curacy or  correctness  of  the  method.  It  will  be  shown  that  under 
certain  conditions  the  average  flow  model  results  are  not  repro- 
ducible due  solely  to  subjective  physical  interpretation  of  the 
lubrication  model. 

REYNOLDS  EQUATION  FOR  ROUGH  SURFACES 

In  this  Section  the  qualitative  appearance  of  a real  rough 

surface  is  described  and  a functional  form  of  the  film  thickness 
resulting  from  significant  surface  roughness  is  hypothesized. 
Statistical  characteristics  relevant  to  the  hydrodynamic  lubrica- 
tion of  rough  surfaces  are  defined.  A general  Reynolds  equation 
for  the  lubrication  of  two  rough  surfaces  is  developed  and  the 
interpretation  of  the  general  equation  for  various  bearing  con- 
figurations is  discussed.  Techniques  for  obtaining  average  values 
of  pressure  and  leakage  from  the  general  equation  are  presented 
and  discussed. 

Film  Thickness  Function 

The  general  appearance  of  a three-dimensional  real  rough 
surface  is  shown  in  Figure  2-1.  The  total  height  H of  any  point 
on  the  surface  is  measured  from  the  arbitrarily  placed  x-y  plane. 
■It  is  hypothesized  that  this  height  can  be  considered  as  being 
composed  of  a nominal  height  h and  a roughness  deviation  5 . Then 
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there  exists  a mean  x-y  plane  precisely  defined  by  h(x,y)  while 


the  deviation  of  the  real  surface  from  this  plane  <5(x,y)  is  a 
statistically  random  process.  In  functional  form: 


H(x,y)  - h(x,y)  + 6 (x,y)  (2-1) 

as  indicated  in  the  figure. 

There  are  several  statistical  quantities  which  characterize 
the  random  roughness  component  5 . By  definition  the  mean  of  6 
is  zero.  An  important  parameter  in  determining  the  magnitude  of 
surface  roughness  effects  on  system  hydrodynamics  is  the  sample 
standard  deviation  a: 


(2-2) 
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Often  the  standard  deviation  is  referred  to  as  the  RMS  roughness 
value.  The  roughness  heights  5 (x,y)  also  have  a frequency  dis- 
tribution f (6)  which  characterizes  the  density  of  points  6 (x,y) 
at  a given  height  value.  This  distribution  function  may  be  in- 
tegrated to  give  the  probability  that  a < 6 < b for  arbitrary 
height  values  a < b: 


b 

P[a  < 5 < b]  - / f(6)d<5  (2-3) 

a 

Aside  from  the  statistics  related  to  the  amplitude  of  5, 

the  surface  roughness  also  exhibits  directional  patterns  in  the 

arrangement  of  surface  peaks  or  asperities.  These  patterns  are 

characterized  by  the  sample  autocorrelation  function  R(X  ,1  ) : 
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where  and  X^  are  called  the  x and  y delay  lengths.  If  6 has 
been  normalized  with  respect  to  a,  i.e.,  if  6 is  replaced  by  o/o 

in  Equation  (2-4) , then  R(X  ,X  ) gives  the  correlation  between 

x y 

the  surface  heights  at  any  two  points  separated  by  X^  in  x and 
X in  y.  For  experimental  and  numerical  purposes,  it  is  easier 
to  use  the  profile  correlation  functions  in  x and  y given  by: 
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Again,  if  6 has  been  normalized  by  o,  then  Rx  and  Rv  give  the 
correlation  between  heights  at  points  separated  by  Xx  along  con- 
stant y or  X^  along  constant  x. 

For  large  delay  lengths  Xx  , the  correlation  is  ex- 
pected to  approach  zero,  while  at  X * 0 the  correlation  becomes 
1.  At  some  intermediate  value  of  X,  the  correlation  drops  to 
0.5.  This  value  is  termed  the  50%  delay  length  XQ  The  50% 
delay  lengths  in  each  direction  are  used  to  define  the  surface 
pattern  parameter  y,  introduced  by  Peklenik  [17]: 


0 . 5x 
^0. 5y 


(2-7) 


This  parameter  may  be  thought  of  as  the  length  to  width  ratio  of 
a representative  asperity  (Figure  2-2) , but  this  is  only  partial 
information  about  the  asperity  shape.  There  are  a variety  of  as- 
perity geometries  for  a given  y such  as  rectangular,  circular, 
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diamond,  etc.  The  asperity  geometry  is  also  a function  of  the  nu- 
merical selection  of  a finite  difference  grid,  and  this  selection 
can  alter  the  hydrodynamic  results  for  the  rough  surface  lubrication. 


It  is  with  these  statistics  in  mind  that  the  general  Reynolds 
equation  for  rough  surfaces  is  derived.  The  resulting  equation  must 
allow  the  two  surfaces  involved  to  have  different  statistics.  The 
solution  for  pressure  and  leakage  from  any  such  equations  must  also 
have  statistical  characteristics  derivable  from  the  surface  statistics. 

General  Reynolds  Equation 

The  development  of  a general  governing  equation  for  pres- 
sure and  leakage  in  a hydrodynamic  system  where  surface  roughness 
is  important  employs  some  simplifying  assumptions.  It  is  usual  to 
assume  laminar  flow  due  to  dimensional  considerations.  Related  to 
this  point  is  the  assumed  validity  of  Newton's  shear  stress  law. 

Body  forces  and  fluid  inertia  are  neglected.  Finally,  the  pressure, 
viscosity  and  density  are  taken  to  be  invariant  throughout  the  film 
thickness. 

Given  the  above  assumptions,  Godet  [15]  has  derived  the 
Reynolds  equation  for  two  rough  surfaces  from  first  principles. 

This  formulation  differs  from  the  classical  approach  because  the 
individual  surface  roughnesses  are  not  lumped  onto  one  surface . 

The  details  of  the  Godet  formulation  were  reviewed  in  the  previous 
annual  report  [1] . 

For  the  geometry  shown  in  Figures  2-3  and  2-4,  expressions 
for  flow  are: 
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Assuming  density  and  viscosity  are 
are  constant.  The  final  result  is 
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Because  of  the  appearance  of  the  quantity  (H^+H^  in  (2-10) , 
Godet  [15]  has  pointed  out  that  the  Reynolds  equation  for  the  rough 
surfaces  agrees  with  the  classical  equation  only  if  = 0.  This 
is  equivalent  to  lumping  the  surface  profile  onto  H2  as  was 
done  by  Christensen,  which  clearly  is  not  correct  if  the  surfaces 
have  different  statistics. 

For  surfaces  having  roughness , it  has  been  hypothesized 
that  the  total  surface  height  H can  be  considered  as  being  com- 
posed of  a nominal  height  h and  a roughness  deviation  5.  Referring 
to  Figure  2-4,  the  surfaces  and  S2  are  now  defined  by: 

H1  * hl  " 51  ' H2  = h2  + 52  (2-11) 


It  is  acceptable  to  lump  the  nominal  film  functions : 

h - h2  - hx  *H  - H2  - H1  = h + S1  + <S2  (2-12) 


The  time  derivatives  of  H^H.,  are  evaluated  with  the  aid  of  Fig- 
ure 2-5.  An  Eulerian  coordinate  transformation  is  assumed  to  be 
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valid  for  the  roughness  functions : 
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(2-13) 


Hence : 
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(2-14) 
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For  the  case  of  interest,  the  point  vertical  velocities 
and  W2  result  from  the  motion  of  the  surfaces  in  the  x- 
direction.  With  the  aid  of  Figure  2-6,  and  W2  can  be 
written: 
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(1  3x  + 3t 

3H2  3H2 

(2-17) 

r2  IT  + It 

The  equations  given  in  (2-11)  to  (2-17)  may  be  used  to  simplify 
the  Reynolds  Equation  (2-10) . without  great  loss  of  generality, 

3 hi 

a steady  nominal  film  thickness  is  assumed  such  that  * 0 and 


3h 

3t 


0.  This  is  the  case  in  most  lubrication  systems.  The  re- 


sult for  Equation  (2-10)  is  then: 


^[<h+w3  H] + 4[(h+w3 !?] " 6n(w  S 


+ 6n (U2-Ux) 


3 <61“<s2> 
3x 


(2-18) 
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Equation  (2-18)  is  the  Reynolds  equation  for  two  rough  surfaces 
as  agreed  upon  by  Godet  [15]  , Elrod  [12]  and  Cheng  [6]  . The  last 
term  in  this  equation  makes  it  clear  that  the  individual  surface 
roughnesses  should  not  be  lumped  together,  as  is  done  in  the 
usual  Reynolds  equation. 

The  application  of  Equation  (2-18)  to  a physical  problem 
should  be  performed  carefully.  This  is  because  the  vertical  veloc- 
ities w 1 and  W2  require  a specific  interpretation  for  each  bearing 
configuration.  A few  examples  of  this  problem  are  discussed  by 
first  simplifying  the  vertical  velocity  expressions  given  by  Equa- 
tions (2-16)  and  (2-17)  with  Equations  (2-11)  to  (2-15)  . With 
3h  3h 

■———  = 0 and  — — = 0,  these  relationships  become: 

3 1 3 1 
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Now  for  a case  of  pure  rolling  where  0 
tive  vertical  velocity  becomes: 


U2  = U,  the  rela- 
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(2-22) 


Using  either  (2-10)  or  (2-18)  , the  right  hand  side  of  the  govern- 
ing equation  becomes: 


3h 

RHS  = 12nu  — 
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(2-23) 


*—■  — * Li. 


27 


This  case  is  compared  to  a situation  of  parallel  plates  where 


0^  » 0^  ■ 0,  but  necessarily 

3hl  _n 
3t  * °1  3x 

at  a2  ax 


0.  In  this  case: 
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0 


(2-24) 
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RHS  - 0. 

In  the  case  of  sliding  cylinders  where  0^-0,  U2 
the  governing  equation  is  solved: 


wL  » 0 
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RHS  - 6n0  ^ + 6nd  
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(2-26) 

(2-27) 
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For  sliding  plates  where  0 ■ 0,  it  is  required  that  WL  - W2 


W. 


0. 


(2-29) 
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which  for  constant  nominal  geometry  implies  that  — 
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(2-30) 
0 . Then 
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Average  Reynolds  Equation 

Equation  (2-18)  governs  the  pressure  of  the  fluid  between 
two  rough  surfaces . Since  the  equation  contains  the  random  sur- 
face deviation  <5,  it  is  not  a useful  expression  unless  £ can  be 
eliminated  in  favor  of  the  statistics  of  6.  If  this  can  be  ac- 
complished,  then  a solution  for  the  expected  value  of  pressure 
can  be  obtained. 

There  are  two  techniques  by  which  results  for  average 
pressure  are  obtained.  In  one  method  the  equation  is  averaged 
prior  to  solution  by  use  of  expectancy  operators.  The  expecta- 
tion expressions  containing  5 are  then  interpreted  in  terms  of 
the  statistics  of  5 and  the  resulting  equation  is  an  average 
Reynolds  equation.  The  solution  of  this  equation  yields  the 
average  pressure  directly.  Alternately  the  governing  equation 
(2-18)  may  be  solved  numerous  times  for  different  but  statisti- 
cally identical  surfaces  and  the  results  for  load  support  and 
leakage  flow  are  averaged.  Both  of  these  techniques  have  been 
used  in  the  past,  the  former  largely  by  Tzeng,  Saibel,  Christensen 
and  Tonder  and  the  latter  by  Patir  and  Cheng. 

Christensen  and  Tonder  assumed  that  surface  S1  has  no  rough- 
ness so  that  £^  * 0.  Then  by  assuming  the  frequency  distribution 
of  62,  Equation  (2-18)  can  be  averaged  directly  by  use  of  expec- 
tancy operators.  The  problem  encountered  in  this  method  is  in 
taking  the  expectancy  of  product  terms  if  the  terms  are  correlated. 

For  example,  considering  the  expected  value  of  the  y-direction 
leakage , 
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(2-32) 
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if  the  terns  (h+6,+5_)J  and  are  correlated.  Christensen  and 
Tonder  found  a way  around  this  difficulty  by  using  simplified 
surfaces  such  as  longitudinal  or  transverse  roughness  patterns. 

Now  in  a preferred  direction  roughness  arrangement  the 
leakage  flow  perpendicular  to  that  direction  does  not  vary  ran- 
domly. This  assumption  allows  the  Reynolds  equation  to  be  aver- 
aged term  by  term  in  a mathematically  valid  fashion.  The  averaged 
equations  for  longitudinal  and  transverse  roughness  given  by 
Christensen  [8,  111  are: 
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(2-34) 


E(p)  is  the  average  pressure,  and  the  expectancy  function  of  a 
quantity  x above  is  defined  in  conjunction  with  the  frequency  dis- 
tribution f of  X! 
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E(X)  - / (x)f(X)dX 


(2-35) 


Given  the  frequency  distribution  f the  expected  pressure  can  be 
found  for  these  roughness  cases.  The  Christensen  approach  to 
three-dimensional  roughness  is  not  accepted  because  both  x and  y 
leakage  flows  are  not  treated  as  random  variables. 

Chong  and  Patir  [6]  have  developed  the  average  flow  model 
which  is  applicable  to  any  generally  non-isotropic  three-dimen- 
sional roughness  structure.  Much  to  their  credit,  the  model  is 
practical  and  the  results  agree  with  the  limiting  cases  given  by 
the  Christensen  Equations  (2-33),  (2-34).  However,  it  is  not 
clear  how  Cheng  and  Patir  derived  their  expected  flow  expression 
for  q^,  and  the  method  will  be  reviewed  here. 

A control  volume  with  base  AxAy  and  height  H * Hj  - is 
selected  to  analyze  the  expected  flow  in  a bearing  segment  (Fig- 
ure 2-7) . It  is  assumed  that  AxAy  is  large  enough  to  contain 
numerous  asperities  but  small  relative  to  the  bearing  dimensions . 
The  unit  flow  in  the  x and  y directions  is  given  by  the  averages: 


y+Ay 


x+Ax 


_ ^ 2 1 - 1 r 

Sc  ” 17  ^ Scdy  ; V* 


(2-36) 


where  q^  and  q^  are  given  by  Equations  (2-8) , (2-9) . Cheng  and 

Patir  define  the  pressure  flow  factors  A , $ and  the  shear  flow 

x y 

factor  * such  that 
ys 
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(2-38) 
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where  the  bars  over  the  symbols  indicate  average  or  expected 
values.  It  should  be  noted  that  EO^-H^  =*  E(H)  f h because  of 

the  possibility  of  contact.  As  presented  by  Patir,  4 and  A 

x y 

compare  the  average  flow  in  a rough  bearing  to  that  of  a smooth 
bearing.  The  term  containing  the  $ factor  does  not  seem  to  re- 
sult from  rigorous  evaluation  of  (2-36) , rather  this  term  is  in- 
cluded to  intuitively  account  for  additional  flow  due  to  sliding. 

Next  a mean  flow  balance  is  performed  on  the  control 
volume  resulting  in  the  continuity  expression: 


35  35  3 (H2-Hl) 

3x  3y  w 3t 

Substituting  the  average  flow  expressions  (2-37),  (2-38)  into 
(2-39) , the  average  Reynolds  equation  is  obtained: 
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(2-40) 


Now  a is  the  standard  deviation  of  the  combined  roughness  5^  + S2: 


2 . 2 
ai  + a2 


(2-41) 


As  h /<3  becomes  large,  the  probability  of  contact  becomes  minimal 
and  surface  roughness  effects  become  small.  This  means  Equation 
(2-40)  should  reduce  to  the  smooth  bearing  equation  which  requires  : 
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As  will  be  shown  later,  the  flow  factors  as  defined  depend  only 
on  the  statistics  of  the  rough  surfaces.  These  factors  are  solved 
by  generating  artificial  surfaces  to  solve  the  general  Reynolds 
equation  and  then  by  averaging  the  results  for  4,  $ and  $ . 

x Y s 

Comparing  Equation  (2-40)  to  the  Christensen  expressions 
for  two-dimensional  roughness,  the  limits  on  4>x  and  for  ex- 
treme surface  conditions  are  predicted: 
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Transverse 
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r -,i-i 

ran 

Longitudinal 

(2-43) 


These  limits  were  verified  by  Cheng  and  Patir  [6 ] , assuming  0^  - U2 . 

Given  the  evaluation  of  $ , <p  and  $ for  a given  nominal 

x y s 

film  thickness  h,  height  distribution  f and  surface  pattern 
parameter  y , then  Equation  (2-40)  governs  the  expected  pressure 
and  is  widely  applicable  to  bearing  problems  with  real  surface 
geometry.  The  advantages  over  the  Christensen  approach  are 
significant,  so  that  the  Christensen  approach  is  not  used  herein. 


Discussion  of  the  Average  Reynolds  Equation 

When  obtaining  the  average  Reynolds  equation  proposed  by 
Patir  and  Cheng  (2-40) , it  was  pointed  out  that  the  term  containing 
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iji  was  not  rigorously  defined.  In  an  attempt  to  define  more 
precisely,  the  assumptions  made  by  Patir  and  Cheng  about  the 
average  flow  Equations  (2-37,  2-38)  may  be  restated.  Writing 


Equation  (2-36)  in  terms  of  expectancy  operators: 
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Comparing  (2-44)  and  (2-37),  the  following  hypothesis  results: 
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(2-45) 


Equation  (2-45)  implicitly  states  that  $x  and  <j>s  are  independent. 
This  assumption  is  checked  later.  In  this  form  it  is  evident  that 
the  term  containing  4>g  is  included  to  relieve  the  dependence  of 
(J>x  on  relative  velocity  O^-l^) . The  problem  of  expanding  the  left 
hand  side  of  (2-45)  has  been  a major  barrier  to  rough  surface 
analysis.  Hence  the  significance  of  this  hypothesis  is  funda- 
mental, and  its  validity  must  be  verified. 

The  derivative  of  $g  can  be  expressed  explicitly  for  the 
case  of  pure  sliding  (U^-U^-  j Ug)  where  the  nominal  pressure 
gradients  are  made  zero  =*  0 j . Using  these  simplifica- 

tions to  write  the  general  Reynolds  equation  in  average  form 
(2-18)  the  following  results: 


3<t>  . 3E(6  -<5,) 

S 1 1 6 

~3x  a 3x 

Alternately  the  average  equation  given  by  Patir 
duced  using  the  same  case : 


(2-46) 

(2-40)  may  be  re- 
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Equations  (2-46)  and  (2-47)  are  mathematically  equivalent  assuming 


toat  ff  * f£  holda* 


Assuming  that  is  independent  as  stated  above,  some  in- 


sight into  the  significance  of  $ in  relation  to  surface  roughness 


may  be  obtained  by  closely  examining  Equation  (2-46)  . It  may  be 
assumed  that  the  roughness  functions  are  truncated  by  asperity 
shear  so  that  E(61)  + 0 and  E(52)  + 0,  assuming  also  that  elastic 
deformation  is  negligible.  If  both  surfaces  are  moving,  the  sur- 


face statistics  are  constant  because  it  is  expected  that  6^  and  <$2' 


will  change  only  very  slowly  in  time  with  wear.  Hence  ECd^)  and 

3<*‘s 

E(52)  are  constant,  meaning  that  • 0.  Then  the  only  term  on 


the  right  hand  side  of  the  average  Reynolds  equation  that  has 


...  3h 

significance  is 


In  the  Patir-Cheng  work,  a different  assumption  appears 
to  have  been  made.  It  appears  that  asperities  are  allowed  to 


spring  back  after  passing  through  contact.  Thus  the  essential 


£ X 

relationship  -r£  - -0  v—  would  not  be  correct  because  of  an  ad- 

d t cX 


ditional  5 (t)  due  to  elastic  deformation.  The  statistics  of  6 


3<p 


and  6 2 also  become  unclear.  For  two  moving  surfaces  -jj-  exist 
rather  than  become  zero  as  in  the  case  above.  Clearly,  more  work 
needs  to  be  done  on  this  aspect  of  the  problem. 

Assuming  the  truncation  model  above , some  of  the  physical 
implications  may  be  illustrated  by  considering  a few  examples. 
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In  the  case  of  two  rotating  cylinders,  one  would  expect  that  the 
roughness  functions  would  take  on  fixed  characteristics  which  would 
change  only  slowly  with  the  progression  of  wear  as  above.  Hence 
the  term  ECS^-^)  is  constant  with  x and  the  average  equation  for 
this  case  becomes : 
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(2-48) 


This  case  is  compared  to  the  case  of  two  parallel  plates  both  of 
which  are  moving.  The  nominal  film  thickness  h is  constant  and 
the  vertical  velocities  W^,  V«2  are  zero.  Then  the  average  equa- 
tion becomes 
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Considering  the  case  of  one  rotating  and  one  fixed  cylinder, 
it  is  physically  reasonable  to  assume  that  the  moving  surface  must 
take  on  fixed  roughness  characteristics.  The  fixed  surface  may 
have  E(<5^)  * f(h),  but  not  a function  of  time,  if  = 0.  Hence 


the  average  equation  becomes: 
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The  basic  debate  which  is  highlighted  in  this  discussion 

concerns  the  physical  existence  of  <J>  . It  has  been  argued  that 
3<j>  3 

for  some  cases  is  zero  and  hence  it  should  not  appear  in  the 
average  Reynolds  equation.  However,  it  does  seem  that  <kg  is  a 
necessary  term  for  evaluating  flow  in  a bearing  segment.  From 
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these  results,  great  care  must  be  used  in  applying  the  average 
flow  model.  Certainly  the  precise  character  of  ^ and  $2  must  be 
understood  before  the  correct  equation  can  be  chosen. 


determination  of  flow  factors  by  simulation 

In  this  section  the  nrcdels  used  by  Patir  end  Cheng  (6.18! 

to  solve  the  presents  flow  factors  ♦ *nd  the  shew  flow  factor 

♦ a are  renewed.  A procedure  to  test  the  ind.pMd.nc.  of  and 

/ i5  detailed.  A Mdel  for  testing  the  bounds  independence  of 

/and  * is  also  proposed.  Two  sethods  to  nueerioally  sinolat. 

x y 

the  flow  models  are  discussed  and  compared. 


Model  to  Solve  Flow  Factors 

A bearing  segment  having  dimensions  Lx  * is  selected 
which  is  small  relative  to  the  bearing  dimensions,  yet  large 
enough  to  contain  numerous  asperities  (Figure  2-8)  . As  will  be 
described  later,  artificial  surface  roughness  functions  $L  and  S2 
are  generated  with  known  frequency  distribution  and  surface  pattern 
parameter  Y to  describe  the  total  film  thickness  in  the  bearing 
segment  at  each  point.  The  segment  has  a constant  nominal  film 
thickness  h - h2  - hr  Hence,  the  governing  equation  for  pressure 
in  this  segment  from  Equation  (2-18)  appears  as: 
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This  equation  can  be  formulated 


numerically  given  the  selection  of 


a finite  difference  network  as  will  be  shown  later.  One  boundary 
condition  is  implicit  in  this  formulation,  that  being  that  no  flow 
is  allowed  at  contact  points.  Hence,  the  film  thickness  expression 
(h  + 5 + $2)  is  set  to  zero  if  this  quantity  should  artificially 

be  generated  negative.  The  remaining  boundary  conditions  are  dic- 
tated according  to  which  flow  factor  is  to  be  found. 

The  pressure  flow  factor  <fr  is  found  by  choosing  boundary 
conditions  which  eliminate  4 and  4 from  the  average  leakage 

y 3 

expressions  (2-37),  (2-38).  Shear  flow  factor  $g  is  eliminated 

by  considering  the  case  of  pure  rolling  * °2  * Factor  ^y 

35 

is  eliminated  by  setting  the  nominal  pressure  gradient  * 0. 

Hence,  the  governing  equation  and  boundary  conditions  for  the  $x 
problem  become: 

■h  [<»*w3  HI  * H [<h*vs2>3  HI  - 0 • <2-52> 
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(2-53) 


Given  a solution  to  the  above  equation  with  boundary  conditions 


for  a given  surface,  is  found  by  eliminating  qx  from  Equations 
(2-8) , (2-37)  and  (2-38)  . Th«  result  for  4>x: 
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where : 


Pb  " Pa 


(2-55) 


The  numerator  of  Equation  (2-54)  is  evaluated  from  the  solution  to 
(2-52)  and  only  one  line  integral  is  necessary  to  find  E(qx)  since 
y flow  is  suppressed  at  the  boundaries.  It  should  be  noted  that 


the  gradient  could  also  be  made  zero  by  using  p(x,0)  ■ p(x,L^)  =■ 

p . However,  Patir  did  not  use  this  condition  and  the  difference 
c 

between  these  boundary  conditions  indicates  that  $ and  4>  may  have 

x y 

a boundary  dependence.  This  can  be  checked  by  a model  to  be  des- 
cribed later  in  this  chapter. 

Since  4>x  is  a random  variable  as  a consequence  of  the  rough- 
ness functions,  it  is  necessary  to  solve  Equations  (2-52) -(2-54)  sev 
eral  times  using  different,  but  statistically  identical  surfaces. 

The  value  of  <(>x  is  then  taken  to  be  E($x),  and  statistical  inference 

concerning  the  range  of  <p  can  be  made.  It  is  observed  that  <p  is 

x x 

independent  of  velocity  U since  this  factor  does  not  appear  in  the 
formulation,  and  that  the  value  of  the  nominal  pressure  gradient 

3p 

is  also  not  a factor  since  the  pressure  is  effectively  normalized 
by  this  quantity  in  the  expression  for  $x-  Hence,  $x  is  a function 
only  of  the  nominal  film  thickness  h and  the  roughness  statistics. 
Also,  since  the  quantity  (<5^  - 5^)  is  eliminated  from  the  govern- 
ing equation,  only  one  surface  is  required  so  that  $x  * $x(h/a,Y) . 
Finally,  it  is  noted  that  since  the  method  of  solving  would  be 

identical  to  that  for  $x>  it  is  unnecessary  to  solve  The  dif- 

ference here  is  that  y-flow  sees  a different  surface  pattern  than 


$ (h/o i y ) - ♦ (h/o,l/y) 
y * 


(2-56) 


The  shear  flow  factor  *g  is  solved  by  eliminating  *x  »d 

* from  (2-40) . The  *g  term  is  retained  by  selecting  ^ / 02  «d 

J , .0  . i i.  selected  for  convenience.  This  problem  becomes: 
U1  2 2 

± (h.5^2>3||  ♦ £ Ji 


a(<5l“<52> 


. (2-57) 


(a)  p(0,y)  - p(ix-y)  * Pa 

0 


“»•  $ 


y-0,L 


(2-58) 


(2-59) 


and  there  is  no  flow  allowed  at  contact  points.  Solving  for  ♦, 
from  (2-8)',  (2-37)  and:  (2-38)  : 

/•Lx  rLy  -(h+6.+6_)  so 
2 I / . 1 .2—  j£  dv  dx  . 

*s  • a CTL "T  J Ja  3x 

s x y o o 

lqain.  ♦ , wat  M avarapad  ovar  savaral  solution,  to  obtain 
It  is  MW  from  (2-57)  that  tha  prossura  is  proportional  to  n 0,. 
and  hano.,  ♦ , dapands  on  h.  7 7,.  d,.  and  o,  as.win,  6,  and  5, 

have  the  same  frequency  distributions. 

Th.  indapondanea  of  and  ♦,  «Y  b.  ohwkad  by  alienating 
only  * fro.  tha  .isolation.  Tha  potarning  aquation  b.co~s  = 


a(«,-«2) 


no.-u 

& HI  * &Kw3  *1  ■ 6n,Dj",i>  “ 57 


(2-60) 


subject  to  boundary  conditions : 
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12-6J 


(b)  p(Lx,y)  = pb 
(C)  3y  y=0,Ly  _ ° 


Elimination  of  qx  from  (2-8),  (2-37)  and  (2-38)  now  becomes: 


jA  - (± 1 IeV  , -=i-  P {h+&L62^  3o 

2 yi2n  3xy^x  LxLy  ,/o  i2n  3x 


dy  dx  . (2-6] 


The  solution  of  (2-60) -(2-62)  will  yield  a line  with  slope  f given 


h 3£ 

6n  U a 3x  ' 
s 


The  hypothesis  to  be  tested  is  that  if  <j>x  and  $g  are  independent 
then  all  solutions  of  (2-62)  for  any  f intersect  in  a common  point 
This  point  should  be  the  solutions  for  <J>x  and  d>g  predicted  by  the 
previous  models.  Since  the  quantities  $x  and  are  random,  the 
right-hand  side  of  (2-62)  actually  has  upper  and  lower  limits 
given  by  the  confidence  interval  on  this  quantity  derived  from 
its  mean  and  standard  deviation.  Hence,  the  lines  given  by  (2-62) 
are  averaged  over  several  solutions  and  actually  plot  as  regions 
bounded  by  two  parallel  lines.  Hence,  a more  precise  hypothesis 
of  independence  requires  that  there  exist  a region  in  the  $x~$g 
space  that  is  a subset  of  all  regions  of  solution  to  (2-62) . 


■k  (<h+ii'Hi2|3  HI  * H ['h*sirt2)3  H 


(2-64) 


Eliminating  q^  and  q^ 


(a)  p(0,0)  ■ p^' 

(b)  p(0,I.y)  - p2 

(c)  P t^jc'O)  “ p4 

(d)  p(Lx,Ly)  - P3 

from  Equations  (2-8) , (2-37)  and  (2-38) ; 
y (h+<S.+<5,) 


(2-65) 


m 11  ^ 1 


1 2 
12 


dy  dx 
3x  1 


-h3  3p 
12  3x 


(2-66) 


where : 


if1  r 

A.  A X 


L*  <h*W  3 


12 


3y 


dx  dy 


3d 

3x 


-h3  3p 
12  3y 


93+P4'V£2  3d  . -VP3"PrP4 


2 L 


3y 


2 L 


(2-67) 


(2-68) 


The  relative  magnitudes  of  and  are  not  a factor  in  the  solu- 
tion of  (2-64) -(2-66)  for  and  $y.  This  is  because  for  a given 


nominal 


|2.,  is  independent  of  because  it  is  a linear  effect. 


The  same  is  true  of 


given 


Ty  ' 3x" 

to  check  the  effect  of  boundary  conditions  on 

Equation  (2-56) . 


However,  this  model  can  be  used 
and  also  to  test 
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1 


Numerical  Formulation  of  Flow  Models 

All  of  the  models  mentioned  require  artificial  surface  gen- 
erations of  5^  and  with  a given  frequency  distribution  and  sur- 
face pattern  parameter  and  mean  zero.  Since  many  engineering  sur- 
faces are  approximately  Gaussian  distributed,  only  these  surfaces 
will  be  considered  herein.  However,  the  phenomenon  of  asperity 


shear  can  lead  to  a truncated  distribution  at  a significant  height 
level  and  this  would  change  the  flow  characteristics  significantly. 


Hence,  the  Gaussian-truncated  distributions  will  also  be  considered 
(Figure  2-9) . 

Artificial  surface  generation  begins  with  the  generation  of 

an  array  2.  . for  i * l,...,n  and  j * l,...,n  such  that  2.  .is 

i.J  x y i,3 

normally  distributed  with  me  am  zero  and  standard  deviation  a.  The 
desired  surface  pattern  parameter  is  imposed  on  the  array  by  a 
transformation  provided  by  Patir  [IS]. 


^ ui 

5i'j  = ™ k?i  £ Zi+k'j+*  ' 


(2-69) 


The  resulting  array  5^  ^ has  linear  profile  correlation  functions 

2 

R and  R such  that  R (0)  * R (0)  = a , R (nAx)  =*  0,  R (mAy)  * 0 
x y x y x y 

where  Ax  and  Ay  are  the  grid  spacings  over  which  the  array  5^  ^ is 
superin^osed.  Since  the  functions  R^  and  R^  are  linear,  the  50% 
delay  lengths  are  given  by: 


i hAx  , _ mAy 

5 x “ 2 r ^ q 5 y 2 


(2-70) 


Hence,  the  surface  pattern  parameter  becomes: 
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n Ax 

m &y 


The  parameter  y is  then  controlled  by  selection  of  n and  m. 

Having  generated  5,  and  5 discreetly  over  the  grid  N x n 
x 4 x y 

it  is  necessary  to  form  the  film  thickness  array.  From  Equation 


(2-11)  : 


h + <S^  + 

i. j ifj 


For  the  Gaussian- symmetric  distribution  the  contact  condition 

iaplies  that  if  H.  . < 0,  set  H.  . = H . where  H . is  some  mini- 
i , j i , j min  min 

mum  possible  film  thickness  greater  than  zero.  The  probability  of 
contact  can  be  found  using  the  Gaussian  density  function  f^($) : 


P[H,  . < 0] 
i » j 


f A < £ ■ 

go  o 


The  truncated  distribution  has  a slightly  different  contact  condi- 
tion. First,  a level  of  truncation  rh  is  selected  (Figure  (2-10) 
Based  on  the  directions  in  which  6^  and  $2  are  defined  in  the 
sketch,  <S^  and  $2  are  first  truncated  according  to: 


S < ~rh  set  = ~rh  + — H ? 

li,j  h.j  2 mn 


If  62  < -(l-r)h  set  $2 

i»  j i » j 


- (l-r)h  + ~ H . 

2 min 


H^  is  then  formed  according  to  (2-72)  and  the  minimum  value  of 

H.  . will  be  H . . Now  the  probability  of  contact  is  the  product 
i , j nun 

of  two  probabilities: 


It  should  be  noted  that  a value  of  r = 1 reduces  the  area  of  contact 
bv  Pi  < 0 or  0.5  of  the  value  for  the  Gauss ian-symme trie  dis- 

L 2i, j J 

tribution.  Hence,  the  truncated  distribution  will  be  characterized 
by  less  contact  than  the  symmetric  distribution.  It  is  also  true 
that  for  a given  h/a,  the  value  r = 0.5  produces  the  most  contact 
and  contact  decreases  with  increasing  or  decreasing  r from  0.5. 

This  situation  of  less  contact  in  the  truncated  case  is  illustrated 
in  Figure  2-10.  The  dotted  lines  which  have  been  truncated  would 
have  otherwise  produced  contact  in  some  cases  for  a symmetric  dis- 
tribution. 

The  next  step  in  numerical  formulation  of  the  model  involves 
the  selection  of  a grid  network  to  serve  as  a basis  for  finite  dif- 
ference equations.  Two  such  networks  are  shown  in  Figures  2-11  and 
2-12 . Figure  2-11  uses  diagonal  elements  which  are  believed  to  be 

the  ones  selected  by  Patir  [18]  for  his  work.  The  next  figure 
depicts  rectangular  elements  which  form  the  basis  for  the  work  herein. 
Each  of  these  element  selections  have  a valid  physical  interpreta- 
tion and  will  be  discussed  separately. 

Patir' s method  places  the  pressures  in  a rectangular  grid 
P.  , for  i = 1,...,N  - 1 and  j = 2 ,..,N  - 1.  The  film  thickness 
function  is  generated  at  both  the  pressure  or  grid  points  and  at 


© 


Bounds  pressure  control  volume 

Bounds  region  of  constant  film  thickness  H 
Pressure  node 


Figure  2-12.  Grid  Network  for  Corner  Point  Method 


the  half-grid  points  resulting  in  H.  . for  i - 1,...,2  N - 1 

Z.  / J X 

and  j * 1 , . . . ,2  N -1.  Next,  two  separate  arrays  HX.  . and 
y i *3 

HY.  . are  selected  from  H.  , at  the  half  grid  points  only: 
i , 3 

HXi,j  " H2i-l,2j  ' 1-1 Nx  j-l,...,Ny-l 

(2-76) 

“i.j  " H2i,2j-l  ' 1 " 1 Nx~  1 j " 1 Ny  ‘ 

H.  .is  generated  with  X » nAx  and  X ■ mAy  so  that 

i,j  0.5x  0.5  y 

nAx  mAy 

HX.  . and  HY.  . have  the  desired  X-  —r—,  X _ ■ • 

1,3  1,3  0*5  x 2 0.5  y 2 

In  order  for  this  element  to  work  the  continuity  equation 
written  for  a control  volume  at  a pressure  node  in  finite  differ- 
ence form  must  correspond  to  the  governing  equation  for  the  model 
under  consideration.  Continuity  for  the  control  volume  may  be 


written: 


, Jh.  JH 
3x  3y  " ' 3t  * 


(2-77) 


It  is  only  convenient  to  analyze  for  the  Patir  element  under 
certain  circumstances.  For  the  case  of  pure  rolling  where 
0 • U2  * 0,  select  At  such  that: 


Ax  - UAt  . 


(2-78) 


Then  t—  is  solved  by  a central  difference  in  time: 

H - ^ + J “I  - - i it)] 

* w [B(x  ’ i - Six  . j ix.y.t)]  . (2-79) 

Since  both  surfaces  are  moving  in  the  same  direction,  the 


expression  in  (2-79)  becomes 


i.j  " Ax  (HXi,j  " HXi+l,j)  ' (2-80) 

When  the  remainder  of  Equation  (2-77)  is  finite  differenced,  it 
is  found  that  the  velocity  flow  terms  on  the  left  side  cancel  the 
■3^  terms.  Hence,  the  right  side  becomes  zero  so  this  element 
correctly  models  the  governing  equation  for  pure  rolling  given 
in  (2-64)  . 

For  pure  sliding  where  U --0  «-Jd  , is  evaluated 

1 i 2 S dt 

as  follows : 

Ax  » j 0s  At  . (2-81) 

ft  " ^ [H(t  + T,X'y>  “ H(t  * 2 At'x'y)j 
0 

* 2/Sc  II  h + 6i(t'x  + J Ax'y)  + 62(t,X  " I Ax'y,j 

- [h  + 51(t,x  - j Ax,y)  + <S2(t,x  + y Ax,Ay)jl.  (2-82) 


3H 

3t 


Equation  (2-82)  assumes  surface  2 moving  in  the  positive  x- direction 
as  implied  by  (2-81).  Converting  to  indices: 


£-*[(» 


+ 5 


1i+l, j 


+ <5. 


>( 


h + 5 


1. 

i.l 


+ 5 


2i+l 


,)!• 


(2-83) 


Although  the  h terms  in  (2-83)  could  be  cancelled,  it  is  easier 
to  enforce  the  contact  conditions  if  the  expressions  in  parenthesis 
are  stored  in  separate  arrays  Hx+  and  Hx~.  The  result  of 
Equation  (2-77)  for  pure  sliding  : 
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Equation  (2-84)  correctly  represents  the  equation  for  pure  sliding 

given  by  (2-57) . It  is  emphasized  that  for  the  Patir  element 

cannot  be  evaluated  in  the  limit  of  it  + 0 since  At  is  finite. 

Hence,  only  the  configurations  of  pure  rolling  (U^  ■ 0^  ■ &)  80(1 

pure  sliding  (U.  - -0,  - - 4 0 ) can  be  analyzed. 

The  second  grid,  termed  the  corner  point  method,  is  based 

on  an  idea  by  Castelli  [19]  and  is  sketched  in  Figure  2-12.  Again, 

the  pressure  array  is  rectangular  and  described  by  -or 

i - 1,...,H  - 1 and  j - 1,...,K  • 1.  In  this  arrangement  the 

film  thickness  array  is  generated  only  at  the  element  midpoints 

and  becomes  H.  . for  i “ 1,...,N  and  j “ 1,...,N  . 

x r 3 x * 

3H  • 

The  analysis  of  -r—  for  the  control  volume  shown  in 

Figure  2-13  is  single  and  general.  With  the  aid  of  Figure  2-13, 

the  volume  V of  the  control  volume  at  time  t is  given  by: 


V(t)  - (Hi#j  +H.fj+i 


+ H 


i+1, j+l'  4 


Mow  subjecting  each  surface  to  a small  time  increment  At  and 
assuming  U > 0. , the  volume  at  time  t +■  At  can  be  written: 


V(t.At)  - j+1l  f <f  * V*> 

* f ^ ' V*> 

*1(>vaJ*KwaJ 

if  <w“  • 

J (2-86) 

The  general  expression  for  t~  is  found  by  taking  the  limit s 


av  ..  V(t+At)  - V(t) 

— ■ lim  

3t  At-O 


>1 


fi  (C^)  (HTR.  # j+HTR.  f j+1)  (H.  # >H.  f j+1>  -Oj  (H.+1  § j+1 

(2-87) 


where: 


HTR 


i*  j 


h + 6 


1, 


+ 6. 


i+1,  j 


i»  j 


Then  ~ is  found  from: 
3t 


3h  _ i av 
at  * AxAy  at 


(2-88) 


(2-89) 


Thus,  it  is  necessary  to  form  the  additional  array  HTR^^  for 

i - 1 sx  - 1 and  j - 1,...#N  prior  to  discarding  5X  and  «2. 

The  contact  conditions  for  HTR  are  the  same  for  H. 

Continuity  for  the  control  volume  can  be  written  in  terms 
of  explicit  mass  conservation,  resulting  in  a completely  general 
Reynolds  equation  which  represents  Equation  (2-51).  After  simpli- 
fication this  formulation  becomes: 
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wmmm 


A.  ,P.  . ■ B.  .P.  . . + C.  .P.  . . + D.  .P.  . , + E.  .P.  . , + F. 

i.j  1,3  i,j  1-1,3  1,3  1+1.3  1.3  1.3-1  i.j  1.3+1  l.j 


. ,„3  , .,3  , 3 3 . ,Ay  Ax. 

A.  . * (H . . + H.  , ,+H.  . , , + H . , . , ) (■— - + r—5  , 

1,3  1.3  1+1.3  1.3+1  i+l, 3+1  Ax  Ay 


B.  . * (H3  . + H?  ) & , 
1,3  1,3  1,3+1  Ax 


C. 


i,j  " (Hi+l,j  + Hi+1, j-l)  A*  ' 


D.  - (H3  . + H3  .)  7^ 
i.j  1,3  1+1,3  Ay 


Ei,j  * (Hi,j+l  + Hi+l,j+l>  Ay  ' 


Pi,j  ’ -12  n(°2  - °l)Ay  [HTRij  + HTRi,j+l 


(2- 


It  is  observed  that  Equation  (2-90)  is  completely  general  and 
reduces  to  the  required  equations  for  pure  rolling  and  pure  slid- 
ing. An  additional  advantage  lies  in  the  fact  that  any  combina- 
tion of  rolling  and  sliding  may  be  solved  with  the  same  equation. 

Patir  has  contended  that  the  film  thickness  at  the  half- 
grid points  should  be  defined  directly  rather  than  averaged  because 
averaging  changes  the  statistics  of  the  surface.  While  it  is  cer- 
tainly true  that  the  half-grid  film  heights  would  have  to  be  aver- 
aged if  the  roughness  was  generated  at  the  pressure  points , it  is 
not  necessary  for  the  comer  point  method  just  described.  Granted 
that  the  finite  difference  equation  appears  to  have  average  terms 
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on  the  left-hand  side,  it  is  not  the  case  in  terms  of  pure  physi- 
cal reasoning.  This  method  has  advantages  of  smaller  roughness 
grids  required,  only  a single  array  describes  each  surface  and 
the  method  is  general. 

Aside  from  the  method  of  formulation,  there  are  other  dif- 
ferences between  the  Patir  method  and  the  corner  point  method  (CPM) 
that  have  physical  meaning.  The  Patir  contact  points  are  half  as 
large  as  the  CPM  and  there  are  twice  as  many  contacts  for  a given 
h and  delay  correlation  values  n,m.  For  sufficiently  large  grids 

N x n this  is  not  a critical  factor.  It  is  also  noticed  that  a 
x y 

single  Patir  contact  only  blocks  leakage  in  one  direction  only. 

This  is  not  realistic.  The  CPM  contact  obstructs  leakage  in  both 
directions,  and  it  would  be  reasonable  to  assume  that  the  CPM  would 
have  a higher  resistance  to  flow.  However,  the  path  length  around 
a Patir  contact  is  considerably  longer  than  the  path  length  around 
a CPM  contact  which  offsets  the  higher  flow  resistance  mentioned. 

These  concepts  can  be  illustrated  by  considering  the  example 
flow  problem  shown  in  Figure  2--14.  The  grids  for  the  Patir  method 
and  the  comer  point  method  are  carefully  drawn  so  that  each  bearing 
segment  has  a total  available  flow  area  of  9 AxAy.  The  film  thick- 
ness in  all  elements  is  taken  to  be  unity  except  for  the  cross- 
hatched  regions  which  are  contacts  where  H ■ 0.  Film  thicknesses 
are  indicated  in  the  elements  and  boundary  conditions  are  enforced 

with  p * 1,  p.  » 0.  Leakage  in  the  y-direction  out  the  ends  is 
a o 

prevented  by  placing  contact  points  in  the  flow  path  at  these 
locations.  TVo  contacts  are  located  in  the  Patir  grid  as  shown 


and  one  in  the  CPM  grid  so  that  the  contact  area  in  each  case  is 

1/9.  Solution  of  $ for  the  Patir  method  is  0.778  while  the  CPM 
x 

d>  is  0.846.  This  exercise  indicates  that  the  Patir  method  offers 
x 

a greater  resistance  to  flow  in  a fundamental  sense.  In  part  this 

is  due  to  the  longer  flow  path  associated  with  the  Patir  element. 

For  example,  the  shortest  path  from  P2  to  p^  is  2 Ay  + Ax,  while 

the  CPM  method  allows  leakage  from  P2  to  p4  through  the  path  length 

Ax.  Even  with  the  contact  between  p.  and  p_  in  the  Patir  method  re- 

4 5 

moved,  only  increases  to  0.800,  still  less  than  the  CPM. 

Because  of  the  advantages  of  the  comer  point  method  over  the 
Patir  method  mentioned,  the  CPM  will  be  used  for  the  results  presented 
in  the  next  section.  The  solution  of  the  finite  difference  equations 
is  by  overrelaxation.  Boundary  conditions  are  enforced  explicitly. 

RESULTS  OF  THE  AVERAGE  FLOW  MODEL 

In  this  section  the  pressure  and  shear  flow  factors  for  iso- 
tropic surfaces  possessing  a Gaussian  symnetric  and  truncated  distri- 
bution of  heights  are  presented.  Results  obtained  from  the  corner 
point  and  Patir  method  of  solution  are  compared  to  each  other  and  to 
published  values  obtained  by  Patir  [18]  and  Elrod  [16] . Then  using 
the  comer  point  method  these  parameters  are  presented  for  non-iso- 
tropic surfaces.  The  boundary  dependence  and  independence  of  <p  , 

<(>  and  <p  are  examined.  Flow  factors  for  the  Gaussian  truncated 
y s 

distribution  are  presented.  Finally  some  estimates  of  solution  ac- 
curacy are  made . 


Flow  Factors  for  Gaussian-symroetric  Surfaces 


The  solutions  for  the  pressure  flow  factor  ^ at  y = 1 or 
for  isotropic  surfaces  are  shown  in  Figure  2-15.  Four  results  are 
plotted  on  this  figure,  two  of  which  are  published  solutions  by 
Elrod  [16]  and  Patir  [18] . Solution  points  obtained  from  the  cor- 
ner point  method  (CPM)  and  the  Patir  method  are  bracketed  by  90% 
confidence  intervals  based  on  10  solution  averages.  The  same  set 
of  different  but  statistically  identical  surfaces  was  used  for  each 
value  of  h/o  in  order  to  obtain  smooth  results. 

Comparing  the  Patir  published  result  to  that  of  the  Patir 
method  obtained  herein,  it  is  observed  that  the  published  result 
barely  intersects  the  confidence  interval  of  the  Patir  method 
results.  Several  seed  numbers  generating  different  sets  of  sur- 
faces were  tried  in  an  attempt  to  obtain  better  agreement  but  none 
altered  this  comparison.  It  is  concluded  that  the  set  of  sur faces 
'used  by  Patir  in  obtaining  the  published  curve  may  have  been  some- 
what of  an  improbable  occurrence,  or  that  the  relative  asperity 
size  to  grid  size  used  by  Patir  differed  from  that  used  herein. 

The  latter  point  will  be  discussed  in  some  detail  later.  In  any 
case  the  published  results  are  statistically  possible  and  are 
considered  satisfactory  by  this  investigation. 

The  solution  curve  for  from  the  CPM  is  qualitatively 
similar  to  the  Patir  result  but  higher  as  was  predicted  from  the 
comparative  discussion  on  the  two  methods.  This  indicates  that 
the  assumptions  about  the  shape  and  character  of  a representative 
asperity  are  important  to  the  outcome  of  the  flow  factors.  The 


square  element  chosen  for  the  CPM  offers  less  resistance  to  flow 


than  the  Patir  element. 

It  is  demonstrated  that  the  results  obtained  independently 
by  Elrod  using  a perturbation  technique  [16]  agree  quite  well  with 
the  Patir  results  for  h/a  > 3.  The  method  used  by  Elrod  is  only 
valid  for  large  h/a  and  the  limitation  on  the  minimum  acceptable 
value  is  not  calculated.  Although  the  results  from  the  comer 
point  method  shown  in  Figure  2-15  do  not  agree  with  Elrod,  there 
are  a few  selective  seed  numbers  for  which  agreement  was  achieved 
for  h/a  >_  4*  0116  such  resuit  is  presented  in  Figure  2-16.  Hence 
all  of  the  results  in  Figure  2-15  are  considered  acceptable  for  the 
regime  h/a  >3. 

The  curves  of  Figure  2-15  also  show  an  increase  in  the  90% 
confidence  interval  as  h/o  ->-0.5.  This  behavior  occurs  for  two 
reasons . First,  the  factor  is  normalized  by  the  nominal  leak- 
age (Equation  2-54) . As  h/a  is  reduced,  the  width  of  the  confidence 
interval  increases  accordingly  since  the  denominator  decreases  with 
(h/a)3.  For  example,  the  reduction  in  nominal  leakage  flow  from 
h/a  * 2 to  h/a  = 1 is  a factor  of  8.  Hence  a 2.5%  scatter  at 
h/a  3 2 becomes  a 20%  scatter  at  h/a  = 1.  Second,  actual  flows 
become  small  as  the  nominal  leakage  becomes  small  resulting  in 
some  machine  roundoff  error.  This  is  a minor  source  of  scatter. 

The  non-isotropic  results  for  obtained  by  the  CPM  are 
presented  in  Figure  2-17.  These  results  are  qualitatively  com- 
parable to  the  Patir  results  (not  shown)  and  in  all  cases  the  es- 
timation of  roughness  effects  by  the  CPM  indicate  less  flow 
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-16.  ♦ for  Isotropic,  Gaussian-symmetric  Surface  Using  CPM 


resistance  than  the  Patir  model  predicts.  It  is  also  evident 
that  the  requirement  $ -*■  1 as  h/o  -►  ® is  satisfied  for  all  values 

of  y. 

Finally,  it  is  observed  that  becomes  large  as  h/o 
tends  to  zero  for  all  values  of  y.  This  behavior  was  not  evident 
from  the  published  Patir  results.  In  fact,  must  tend  to 
infinity  as  h/o  + 0 by  definition  in  Equation  (2-54)  and  this 
point  preserves  the  continuity  of  $x  for  the  complete  range  of  y. 

Before  discussing  the  results  for  shear  flow  factor,  some 

clarification  of  the  parameters  affecting  this  factor  must  be 

made.  It  was  mentioned  in  the  previous  chapter  that  4>s  depends 

on  the  nominal  geometry  h/o  and  the  individual  statistics  of  <5  , 

62.  Now  intuitively  #g  is  flow  trapped  in  the  surface  valleys 

2 2 

under  conditions  of  sliding  such  as  0^  ■ -U2 . Hence  if  = o2 
it  is  expected  that  trapped  flow  in  one  surface  will  cancel  out 
with  trapped  flow  in  the  opposite  surface  when  0^  * The  de- 

pendence of  $g  on  and  <j2  is  postulated  and  proven  to  be  [18] : 


♦s  m(-f)  *(h/c,Y2)  "(y)  *(l»/o,Y1>  (2-91) 

where  $ depends  only  on  h/o  and  y • $3  is  solved  by  letting 

Yx  • y2  • y* 

2 . 
a a 

*s(h/<J,y)  - -5 2 (2-92) 


In  results  not  included  herein  the  non-dependence  of  <fg  on  ff^,a2 
is  verified  and  the  dependence  of  $g  on  o^,o 2 is  removed  by  the 
hypothesis  of  $g  in  Equation  (2-91)  . 


w 


2 2 

The  results  for  $ using  isotropic  surfaces  and  a ^ a are 
s 12 

presented  in  Figure  2-18.  It  is  observed  that  the  published  curve 
has  a slightly  different  shape  from  the  Patir  method  curve.  Patir’s 
published  curve  intersects  the  10  point  90%  confidence  interval 
everywhere  except  the  region  from  3 <_  h/o  <_  5 . This  behavior  can- 
not be  explained  without  more  intimate  knowledge  of  Patir' s work. 

Again  the  CPM  solution  behaves  as  anticipated.  Since  this 
method  offers  less  resistance  to  pressure  forced  flow,  it  also 
drags  less  self-induced  flow  in  the  absence  of  a nominal  pressure 
gradient.  The  results  are  everywhere  less  than  the  Patir  result 
and  the  closeness  to  the  Patir  result  is  about  the  same  as  for  4>x . 

For  this  flow  factor  it  is  observed  that  the  scatter  in 
results  does  not  increase  as  h/o  -*■  0.  This  is  because  shear  flow 
is  normalized  in  terms  of  the  sliding  velocity  us  = U2  ~ U1 
(Equation  2-59)  which  is  invariant  of  h/o  . The  size  of  the  con- 
fidence interval  is  related  to  the  magnitude  of  the  flow  factor 
only,  but  is  subject  to  some  roundoff  as  h/o  -*■  0. 

Shear  flow  has  a maximum  value  for  h/o  near  1.25.  This  is 

because  as  h/o  is  decreased  from  6,  the  pressure  gradients,  which 

are  nominally  zero,  increase  due  to  the  surface  roughness  effects. 

However,  when  h/o  becomes  small  enough,  the  increase  in  pressure 

gradient  is  offset  due  to  the  decrease  in  "valley"  volume.  Hence 

as  h/o  -*■  0,  tends  to  a small  value  not  necessarily  0.  As  h/o  -*■  00 

it  is  observed  that  $ tends  to  zero,  as  well  as  3$  /3(h/o).  This 

s s 

is  noted  because  the  actual  quantity  of  interest  is 
3$  /3x  * (3<f>  /3h) (3h/3x)  in  Equation  (2-40) . 

3 a 
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Figure  2-18.  ♦ for  Isotropic,  Gaussian-symmetric  Surfaces 
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$g  for  non-isotropic  surfaces  as  predicted  by  the  CPM  are 
shown  in  Figure  2-19.  These  results  indicate  less  flow  due  to  sliding 
than  the  Patir  results.  The  reason  shear  flow  increases  with  de- 
creasing Y is  that  the  valley  volume  is  greater  for  more  trans- 


verse surfaces  since  the  valleys  are  longer  in  the  direction  per- 


pendicular to  sliding.  For  purely  longitudinal  surfaces  it  is 
expected  that  very  little  flow  would  be  picked  up  since  there  are 
no  transverse  surfaces  to  trap  the  fluid  in  the  valleys. 


Independence  of  Flow  Factors 


As  described  previously,  the  independence  of  (j»x  and  <|>g 
is  tested  by  imposing  both  a sliding  velocity  and  a nominal  pres- 
sure gradient  By  varying  the  relative  magnitude  of  the  sliding 

velocity  and  the  nominal  pressure  gradient,  the  lines  plotted  in 


Figure  2-20  are  determined.  Each  solid  line  represents  the  mean 
f-H3  3p  1 

for  e|'i2—  3^1  311(3  3s  bounded  by  90%  confidence  intervals  based  on 
10  point  averages.  It  is  evident  that  and  $g  are  indeed  indepen 
dent,  for  all  regions  bounded  by  confidence  intervals  intersect  the 
dark  region  which  contains  the  previous  solutions  for  and  $ at 
h/a  * 2,  y * 1. 

The  boundary  dependence  of  $ and  $ is  tested  by  completely 

x y 

specifying  pressure  on  the  boundary,  as  described  in  the  previous 

chapter.  Results  for  y “ 1*  using  Gaussian-symmetric  surfaces, 

are  tabulated  in  Table  2-1.  It  is  first  observed  that  the  values 

<p  and  <j>  do  not  depend  significantly  on  the  ratio  of  the  magnitude 
x y 

of  the  pressure  gradients  applied  to  the  model  as  predicted.  Next, 
it  is  noted  that  ■ <p^  satisfactorily  for  y » 1 as  they  should. 


Table  2-1.  Testing  the  Boundary  Dependence 

of  4 and  $ . y ■ 1. 
x y 


!£  /I £ 

3 x/  3y 

1 

l 

h/a 

2 

3 

00 

$ 

X 

ax 

1.1400 

0.9707 

0.9693 

y 

m 

• 

• 

3 

* 

X 

OB 

1.1381 

0.9650 

0.9656 

*T 

m 

1.0616 

0.9091 

0.9354 

2 

♦, 

m 

1.1372 

0.9621 

0.9638 

y 

s 

1.0683 

0.9262 

0.9464 

1 

X 

m 

1.1138 

0.9536 

0.9583 

1 

* 

y 

ax 

1.0751 

0.9433 

0.9575 

rzar:-- — - -stcc-r 


However,  the  solution  for  does  not  agree  with  that  given  in 
Figure  2-15.  The  flow  resistance  here  is  less  ($x  greater)  than 
previously  predicted  by  a significant  amount.  This  is  because 
when  the  flow  is  obstructed  by  a contact  it  need  not  pass  around 
the  contact  necessarily,  rather  it  may  just  turn.  Hence,  these 
boundary  conditions  provide  a path  of  lower  resistance  to  the  flow 
when  contact  points  sure  encountered  than  in  the  model  used  for  the 
results  of  Figure  2-15. 

The  explanation- just  given  can  be  tested  by  considering  $x 
and  for  non-isotropic  cases.  Two  values  of  y + 1 sure  tabulated 
in  Table  2-2.  It  is  verified  that  $ (h , y)  ■ (h,l/y)  considering 

y * 

confidence  intervals.  Comparing  the  values  in  Table  2-2  to  those 
presented  in  Figure  2-17,  the  flow  is  everywhere  greater  in  the 
completely  pressure  specified  boundary  case.  This  reinforces  the 
argument  that  the  pressure  flow  factors  are  boundary  dependent  and 
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Table  2-2. 


Boundary  Dependence  of  $ 


4>  for 
y 


Non-Isotropic  Surfaces 


1 


Y 

i 

1 

1 

1.5 

h/a 

2 

3 

1/3 

»x 

- 

0.6686 

0.6548 

0.7025 

0.8026 

3 

1.6289 

1.3373 

1.2176 

1.1139 

1 

*x 

38 

1.1138 

- 

0.9536 

0.9583 

3 

1.0751 

— 

0.9433 

0.9575 

3 

1 

*x 

3 

1.6356 

1.3414 

1.2209 

1.1158 

3 

0.8199 

0.7569 

0.7759 

0.8430 

must  be  redetermined  for  each  particular  application.  Although 


it  is  not  investigated  herein,  it  is  reasonable  to  assume  that  $ 
is  also  boundary  dependent  to  an  even  greater  extent  than  4>x  and 
d>  . Since  flow  can  run  down  the  valleys  transversly  in  this 

y 

model,  less  fluid  will  remain  trapped  by  the  transverse  ridges. 


Flow  Factors  for  Gaussian-Truncated  Surfaces 

The  Gaussian- truncated  distribution  of  roughness  heights 
adds  the  additional  variable  r,  the  truncation  level,  to  the 
analysis  of  flow  factors.  This  level  is  variable  because  of  the 
possible  mating  of  a surface  with  higher  shear  strength  to  a sur- 
face having  lower  shear  strength.  The  significant  truncation 
level  of  the  weaker  surface  will  possibly  be  lower  them  the 
stronger  surface  and  the  probable  cross-section  area  at  a given 
height  must  be  considered  to  determine  r for  a pair  of  surfaces. 

Another  aspect  of  truncated  surfaces  is  that  asperity  deforma- 
tion is  no  longer  present.  This  allows  the  flow  factors  to  be  solved 


for  h/cr  < 0.5,  which  was  a limitation  on  the  average  flow  model  pre- 
viously. Absence  of  asperity  deformation  also  means  that  E(S^)  and 
E ( 6 2)  are  time  and  position  invariant,  lending  support  to  the  previous 
arguments  debating  the  average  Reynolds  equation  techniques. 

Before  presenting  the  results  for  flow  factors  from  truncated 
distributions  it  is  pointed  out  that  no  direct  comparison  to  the  sym- 
metric distributions  can  be  made.  This  is  because  even  with  r = 0.0, 
the  lower  surface  (Figure  2-10)  possesses  valleys  which  reduce  the 
contact  area.  If  is  then  made  0,  a truncated  distribution  no  longer 
exists  and  the  situation  resembles  Gaussian-symmetric.  Hence,  the  re- 
sults to  be  presented  stand  alone  and  can  only  be  compared  to  previous 
results  qualitatively. 

The  pressure  flow  factor  <f>x  for  isotropic  truncated  surfaces 
having  variable  r and  * c is  plotted  in  Figure  2-21.  Values  of 
r from  0 to  0.5  are  given  only  since  for  a1  = a $x(h,Y,r)  * 4>x(h,y,l-r)  . 
As  r increases  to  0.5,  flow  decreases  because  the  contact  increases. 
Qualitatively,  the  results  resemble  those  for  symmetric  surfaces,  but 
the  flow  resistance  is  much  less  due  to  less  contact.  It  is  also  true 
that  the  average  value  of  H is  greater  than  the  symmetric  distribution 
because  deformable  asperities  which  would  decrease  H have  been  trun- 
cated. Hiis  explains  why  a y = 1 in  this  case  nearly  corresponds  to  a 
Y * 2 or  3 in  the  previous  case  of  deformable  asperities. 

The  effect  of  y on  with  r * 0.5,  ® o2  is  shown  in  Figure 

2-22.  Again,  the  qualitative  behavior  is  the  same  as  for  symmetric 
surfaces,  but  the  flow  resistance  is  much  lower. 

Results  for  the  $ case  are  much  more  interesting.  The 
2 2 

linearity  of  4>  with  a.  and  a,  proposed  by  Patir  for  deformable 
s 1 6 
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surfaces  is  no  longer  valid.  To  illustrate  this  a graph  of  <j>s 
for  a = a2  and  variable  r is  shown  in  Figure  2-23.  In  the  symme- 
tric density  function  results , <j>  was  expected  to  be  zero  for 
o * a2-  However,  this  is  only  true  for  r = 0.5  in  the  truncated 
distribution.  Other  values  of  r create  valley  depths  on  one  sur- 
face larger  than  those  on  the  opposing  surface  which  produces  a 
net  shear  flow.  The  curves  have  the  same  basic  shape  as  the  pre- 
vious results  for  4>  , but  the  magnitude  of  the  flow  is  much  lower 
s 

for  any  comparable  case.  Figure  2-24  shows  (ji^  for  fixed  r * 0.5 
2 2 2 

and  variable  o . Again,  for  a1  = o2  * 0.5,  no  flow  is  developed. 
But,  as  a * ■*  0.9,  the  valleys  on  &1  become  larger  them  on  &2  and 
flow  picks  up. 

Results  for  using  non-isotropic  surfaces  truncated  at 
r = 0 with  a^2  = cr22  are  plotted  in  Figure  2-25.  This  figure 
yields  no  new  information  since  the  qualitative  trends  are  con- 
sistent with  those  results  of  the  Gaussian  symmetric  distribution. 

Limitations  of  the  Numerical  Formulation 

Several  error  estimates  and  checks  were  made  on  the  numeri- 
cal method.  The  solution  to  the  finite  difference  equations  was 
accomplished  by  successive  overrelaxation.  Hence,  the  convergence 
criteria  involved  a relative  pressure  change  tolerance  t.  Selec- 
tion of  the  appropriate  value  of  t involved  a tradeoff  between 
computer  time  and  accuracy  of  the  solution.  It  was  noted  pre- 
viously that  <p  could  be  solved  at  any  value  of  x,  since  it  should 
be  the  same  everywhere.  Table  2-3  gives  typical  values  of 
<4  at  x - Ax  (I  * 2)  , x = L - Ax(I  » N - 2)  and  * mean  for  a 
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Table  2-3.  Accuracy  Problem  - Necessary 
Tolerance  for  Convergence. 


T 

<p  1 = 2 

X 

♦J1  1 

<j>Jmean 

IBM  360 
CPU 

io“3 

0.5895 

0.9796 

0.8500 

19 

** 

1 

O 

0.7058 

0.7441 

0.7305 

31 

IQ’5 

0.7191 

0.7153 

0.7177 

47 

IQ'6 

0 . 7196 

0.7142 

0.7172 

>56 

specific  surface.  Even  though  the  end  values  fail  to  match, 

-4 

t * 10  yields  less  than  2%  error  in  mean  values.  This  value 
was  used  for  all  results  and  this  error  estimate  was  repeatable. 

A different  sort  of  limitation  was  inposed  by  considera- 
tion of  the  asperity  size.  This  parameter  is  controlled  by  the 
correlation  lengths  n and  m used  to  generate  the  artificial  sur- 
faces (Equation  (2-66)).  Larger  values  of  n and  m create  larger 
contact  points,  but  fewer  of  them  so  that  contact  area  remains 
the  same.  For  a sufficiently  large  grid  the  values  of  n and  m 
should  not  effect  the  flow  significantly.  However,  when  this 
effect  was  studied  for  isotropic  surfaces  (m»n)  the  effect  on 
<Px  was  profound.  Tables  2-4  and  2-5  illustrate  this  study  for 
the  Patir  method  and  the  comer  point  method.  While  there 
appears  to  be  a trend  of  decreasing  as  n(*m)  increases,  this 
was  found  in  some  instances  not  to  be  the  case.  However,  the 
trend  exhibited  by  the  confidence  interval  size  as  it  increases 
with  increasing  n was  found  to  be  always  the  case.  This  clearly 
indicates  that  a large  range  of  results  for  $ is  statistically 
possible  by  varying  the  choice  of  m and  n for  a fixed  grid  size. 
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Table  2-4.  Patir  Method  - Effect  of  Asperity  Sire  on  (Plow 

Factor) , y (Pattern  Parameter)  and  P (Contact  Area) , 

c 


N(-M) 

Grid 

5 

X 

is  26  x 26,  Y ■ 

90%  Cl 

1,  h/c  » 

y 

1. 

90%  Cl 

P 

c 

1 

0.6179 

0.5671-0.6687 

0.9947 

0.9693-1.0200 

16.00 

2 

0.5937 

0.5186-0.6688 

0.9957 

0.9455-1.0458 

15.97 

3 

0.5780 

0.4401-0.7159 

0.9709 

0.8807-1.0611 

15.96 

4 

0.5379 

0.3299-7459 

0.9142 

0.8462-0.9822 

16.15 

Table 

2-5.  Corner  Point  Method  - 
on  $x,  y.  Grid  is  26 

Effect  of  Asperity  Size 

* 26,  y * 1,  h/a  » 1. 

N(-M) 

$ 

X 

90%  Cl 

Y 

90%  Cl 

1 

1.1191 

1.0143-1.2238 

0.9820 

0.9504-1.0135 

2 

0.8696 

0.7786-0.9607 

0.9798 

0.9082-1.0514 

3 

0.7963 

0.6939-0.8986 

0.9772 

0.9021-1.0522 

4 

0.6393 

0.4639-0.8147 

0.9558 

0.8256-1.0861 

Unfortunately , this 

choice  was  limited  by  computer  space  and  so 

the  value 

am*  2 was 

used  throughout, 

with  a grid 

size  of  26  x 26 

This  value  was  selected  because  it  is  felt  that  n*n»l  is  not  a 
well-correlated  surface  and  because  higher  values  of  m and  n 
reduce  the  number  of  correlation  lengths  contained  in  the  grid. 

The  effects  of  varying  the  grid  size  yields  similar  con- 
clusions as  varying  the  asperity  contact  size.  For  a constant 
contact  size  n,m  the  increased  grid  size  effectively  reduces  the 
relative  asperity  contact  size.  An  example  of  this  is  shown  in 
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Grid  Size  Points  Averaged  a 

NXXNY  NAV lx 


10  X 10 

20 

0.5166 

20  x 20 

17 

0.3477 

30  x 30 

17 

0.2837 

o 

X 

O 

14 

0.1418 

The  studies  of  relative  asperity  contact  size  to  grid  size 
indicated  that  larger  grids  with  longer  correlation  lengths  are 
recommended  for  three  reasons.  First,  the  asperity  size  will  not 
influence  the  value  of  ^ if  a sufficiently  large  grid  xs  used. 
Second,  it  has  been  shown  that  the  confidence  interval  on  is 
reduced  for  larger  grids.  Finally,  it  is  felt  that  the  lubrica- 
tion behavior  would  be  more  accurately  modeled  if  several  pressure 
nodes  were  contained  in  a given  contact  instead  of  only  3 or  4 as 
is  the  case  for  m»n=*2.  However,  a network  with  100  x 100  nodes 
and  n*m-10  would  only  be  feasible  on  a larger  computer  and  using 
a different  solution  technique. 
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CONCLUSIONS  ON  ROUGHNESS  MODEL 


The  following  conclusions  have  been  made  based  upon  this  in- 
vestigation. 

1.  The  selection  of  a finite  difference  grid  network  alters 
the  results  from  the  average  flow  model  significantly. 
In®  licit  in  this  selection  is  the  contact  shape  and 
characteristics.  It  is  felt  that  since  the  Patir 
method  element  does  not  block  flow  in  two  directions, 

it  is  not  a good  selection.  Hie  fact  that  a different 
physical  interpretation  of  the  contact  geometry  alters 
the  flow  factor  suggests  that  more  work  is  needed. 

2.  Since  the  corner  point  method  used  herein  intersects 
the  Patir  and  Elrod  results  statistically  for  > 3, 
it  is  concluded  that  this  method  is  valid  in  this  re- 
gime. It  is  also  thought  that  the  corner  point  method 
is  only  one  of  several  possibly  valid  models  of  the 
rough  surface  lubrication.  The  corner  point  method 
offers  a lower  resistance  to  flow  than  the  Patir  meth- 
od so  that  pressure  flow  factors  are  larger.  This 
conclusion  emphasises  the  point  that  the  average  flow 
model  should  not  be  taken  at  face  value,  rather  the 
lubrication  problem  to  be  modeled  should  be  examined 
carefully  first. 

3.  While  it  is  not  evident  that  the  pressure  flow  factors 
are  continuous  in  y from  the  Patir  published  curves, 


84 


this  continuity  is  observed  to  be  necessary  by  defi- 
nition of  $ and  the  results  presented  herein.  Since 

$ becomes  very  large  for  — < 0.5  at  all  values  of  y, 

X 0 

it  is  questionable  whether  or  not  any  results  in  this 
range  are  usable.  Although  it  is  admitted  that  the  de- 
formation model  is  not  accurate  in  this  range,  it  is 
pointed  out  that  the  truncated  model  remains  valid  as 

— -*>  0.  However,  in  either  case  the  question  still 
o 

stands  due  to  large  scatter  of  $ . 

4.  While  the  shear  flow  factor  $ exists  to  preserve  con- 
tinuity in  the  bearing  segment,  it  is  uncertain  whether 

3*s 

or  not  — — should  exist  in  some  cases.  The  arguments 
3x 

3*3 

that  — — = 0 for  a truncated  surface  are  presented 
3x 

in  the  second  section.  However,  later  results  indi- 

es 

cate  that  -rr-  ¥ 0 for  constant  <5..  and  5_  surfaces.  It 
3x1  1 ^ 

is  suggested  that  this  could  be  a boundary  effect. 

5.  The  results  of  Tables  2-1  and  2-2  indicate  that  the 
pressure  flow  factors  are  boundary  dependent.  Sup- 
pression of  the  nominal  y-direction  leakage  rarely 
occurs  as  in  the  model  used  by  Patir.  Typically  there 
is  a zero  nominal  pressure  gradient  applied  across  the 
bearing  width.  Results  from  this  type  of  boundary  con- 
dition indicate  higher  pressure  flow  factors  than  the 
Patir  model  due  to  lower  flow  resistance. 
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It  is  concluded  from  Figure  2-20  that  the  implied  in- 


dependence of  and  is  valid.  Hence  it  is  ascer- 


tained that  both  the  <4  and  $ terms  are  reouired  to 

x s 


expand  E ^-H3  |^.  The  expansion  of  this  term  when  H' 


and  are  correlated  has  been  debated  for  years , and 
3x 


the  hypothesis  proposed  by  Cheng  and  Patir  for  de- 
formable surfaces  is  valid  at  least  for  an  element  with 
zero  flow  at  the  edges. 

In  both  the  Patir  method  and  the  corner  point  method 
it  is  found  that  the  relative  contact  size  to  bearing 
dimension  alters  the  quantitative  flow  factor  results. 
It  is  concluded  that  the  results  presented  by  Patir 
and  herein  should  be  checked  using  larger  grid  sizes 
prior  to  application. 

It  is  felt  that  the  truncated  model  ~ s more  realistic 
than  the  deformable  rough  surface  concept,  at  luast 
for  cases  when  elastic  deformation  is  negligible.  This 
model  necessarily  requires  the  additional  parameter 
specification  r,  or  the  truncation  level.  It  is  also 
evident  from  the  results  that  both  pressure  and  shear 
flow  factors  depend  on  individual  surface  statistics 
for  truncated  distributions.  It  is  concluded  from 


these  results  that  truncated  surfaces  are  characterized 
by  less  resistance  to  flow  due  in  part  to  less  contact 
and  also  to  the  fact  that  E(H3)  is  greater  for  trun- 


cated distributions 


SIGNIFICANCE  FOR  ROUGH  WAVY  SEAL  LUBRICATION 

The  model  developed  in  the  previous  annual  report  (1]  was 
based  on  the  Christensen  [20]  and  Tonder  [21]  result  for  isotropic 
surfaces  where  the  moving  surface  was  assumed  to  be  perfectly 

smooth.  In  essence,  Christensen  and  Tonder  indicate  that 

3 3 3 3 

<j>x  = E(H  )/h  for  isotropic  roughness.  However,  E(H  )/h  >1. 

Comparing  this  result  to  the  curves  in  Figure  2-15  for  isotropic 

roughness  shows  that  the  <f>x  used  in  the  previous  work  is  too  large . 

Flow  resistance  in  the  previous  work  is  larger  than  predicted  by 

the  present  work;  therefore  the  previous  model  underestimates  the 

hydrodynamic  load  support. 

Before  the  amount  of  the  underestimate  can  be  calculated 
based  on  the  new  model,  certain  other  features  of  the  model  require 
evaluation.  In  the  Patir -Cheng  [6]  approach  the  much  debated  but 
real  factor  <j>g  has  been  introduced.  Now  <f>s  is  significant  only  if 
it  changes  with  respect  to  x,  according  to  Equation  (2-40) . Also, 
for  equal  roughness  on  both  surfaces  <j>s  = 0.  For  the  sake  of  avoid- 
ing the  problem  of  interpreting  4>s,  let  it  be  assumed  that  both 
surfaces  have  equal  roughness  or  that  the  roughness  is  constant 
with  respect  to  x.  It  is  quite  reasonable  to  assume  that  rough- 
ness is  constant  on  the  moving  face  and  has  a truncated  Gaussian 
character  due  to  wear.  On  the  fixed  wavy  face  of  the  seal,  it  is 
likely  that  there  is  some  variation  of  roughness  with  x,  but  this 
can  be  assumed  to  be  negligible ,and  it  can  be  assumed  that  the 
roughness  is  truncated.  Finally,  inspection  of  Figure  2-15  for 
Gaussian  surfaces  and  Figure  2-21  for  truncated  Gaussian  surfaces 
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shows  that  ^ = 1 is  a good  assumption  for  equal  roughness  sur- 
faces with  h/a  > 2. 

From  the  above  discussion  it  becomes  apparent  that  under 
certain  conditions  the  effect  of  surface  roughness  on  lubrication 
becomes  insignificant,  that  is,  * 1,  =*  0.  A valid  model  can 

be  constructed  without  having  to  answer  the  many  open  questions 
raised  in  the  analysis  of  the  Patir  Cheng  [6]  method.  Osing  the 
above  simplifying  assumptions , the  wavy  rough  seal  configuration 
must  appear  as  shown  in  Figure  2-26.  In  order  to  have  a finite 
region  of  mechanical  contact  (when  needed  for  load  support) , it 
is  necessary  that  the  wave  be  flattened  as  shown.  This  is  required 
if  a constant  truncation  level  distribution  is  used.  In  reality 
there  may  be  some  variation  in  truncation  level  near  the  contact 
region  particularly  if  the  wave  is  fixed.  However,  for  a traveling 
wave  the  above  assumption  appears  reasonable.  Figure  2-26  also 
shows  the  assumed  distributions  for  the  two  surfaces.  Contact  oc- 
curs when  h * 2 tr-  The  mean  value  of  the  untruncated  Gaussian 
distribution  is  used  to  define  h.  The  truncation  level,  tr  as 
well  as  and  a2>  the  standard  deviations  of  the  roughness  dis- 
tributions are  variables  that  must  be  chosen. 

For  a cosine  wave  function,  h is  given  by 

h =*  h + h cos  n0  (2-93) 

z o n 

h * h when  h > 2 t (2-94) 

Z Z “ " * 

h * 2 tr  when  h < 2 t 
z r 


38 


~~ 


(2-95) 


With  the  above  assumptions  and  converting  to  polar  coordi- 
nates, the  dimensionless  Reynolds  equation  (2-40)  becomes 


- 38  \ 38  M \ c)  38 


(2-96) 


where  h is  the  nominal  film  thickness  as  defined  and  p is  the 
a o 

average  pressure.  Based  upon  Equations  (2-37)  and  (2-38)  and  let- 


ting $ =0  and  <p 

S X 


<py  * 1,  flow  equations  are 


3 _ 

^9<j  “ IF  + 2 [ + E(ha  + ^2a]  ] 

12r 


(2-97) 


3 _ 

-E  3p 
a a 


ra 


12 


(2-98) 


3r 


The  above  dimensionless  quantities  are  the  same  as  in  Reference 
[1] , but  have  been  normalized  using  a instead  of  c.  Thus,  using 
these  simplifying  assumptions,  the  only  term  that  reflects  the 
presence  of  roughness  is  E(<S^  + <S2)  . cav^^  boundary  con- 

dition, Equation  (2-37)  in  Reference  [1],  now  becomes 
.3 


i£! 

3h 


<D 


l©c  = 6r  cosa©  * 


(ha  + E(5la)  ♦ E(S2a))  - (ha  + E(6la)  ♦ S(62o) 


y 


(2-99) 

However,  for  the  given  assumption  that  roughness  is  constant  with 

position,  the  E(5,  ) + E(6.  ) terms  on  the  right  hand  side  cancel 
LO  60 


and  vanish  completely. 
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The  previous  four  equations  along  with  the  straightforward 
boundary  conditions  outlined  in  Reference  [1]  form  a complete  set 
of  equations  for  the  solution  for  pressure  distribution  for  rough 
wavy  seal  faces  under  the  roughness  assumption  outlined.  These 
equations  show  that  for  this  particular  set  of  roughness  assump- 
tions, the  pressure  distribution  is  not  affected  by  the  roughness 
itself.  However,  mechanical  load  support  and  fluid  and  mechanical 
friction  are  affected  by  surface  roughness,  and  a method  of  com- 
puting these  factors  must  be  established  before  a solution  to  the 
wavy  seal  problem  can  be  obtained. 

Mechanical  load  support  is  given  by  the  probability  of  two 
asperities  touching  times  the  total  face  area  times  the  compressive 
strength,  similar  to  reference  [1] . 

p = p P f ( 5 < t ) and  (5  < tr)  ] (2-100) 

m m i i 

a 

Assuming  that  the  roughness  on  the  two  surfaces  is  uncorrelated, 
then 

P C (<S2  < tr)  and  (a  < tr)  ] = P(62  < tr)  x P($L  < t.)  (2-101) 

where 

-t 

r 

P(5  < tr)  = / f<6)d«  (2-102) 

—00 

where  f ( 5 ) is  the  Gaussian  distribution.  P(<52  < tr)  is  the  same 
as  PfS^  < tr)  for  equal  roughness. 

Mechanical  friction  shear  stress  is  given  by  a similar 
relationship. 
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PsP^52  < and  < 


(2-103) 


P„ 


Both 


p and  p are  zero  when  h > 2t  , i.e.,  when  the  seal  does 
m s r 


a a 

not  contact. 

To  calculate  fluid  film  friction  as  in  Reference  [1]  one 

needs  to  evaluate  E A . For  uncorrelated  roughness,  this  is 

H 

given  by 


- izxhr)-// 

mmOO  -.00 

Ohder  the  condition  where  h * 2t  , then  can  become  large,  so  that 

i.  n 

the  minimum  value  of  H is  set  at  A as  in  the  previous  report  [1]. 
Then  the  integral  appears  as  follows: 


00  00 

/,  s«7r-fl5i,5(VdSidS2  (2-105> 

_tr+  2 *tr+2 

For  this  case  of  contact,  f(6^)  and  ffSj)  are  the  Gaussian  distri- 
bution functions,  and  the  probability  of  the  combined  inpulses  shown 
in  Figure  2-26  is  ignored  since  this  is  accounted  for  by  mechanical 
friction. 

When  h > 2tr  + A,  then  the  integral  becomes 


E‘if’  ■ / / (J-1061 

-t  -t  12 

and  the  impulses  on  fCS^)  and  f(62)  must  be  included  since  there 
is  no  mechanical  friction. 

Combining  the  above  relationships  with  the  relationships  and 
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methods  outlined  in  the  previous  report  [1] , one  can  solve  for 


the  pressure  distribution,  friction,  and  leakage  for  this  new 
problem  in  a manner  similar  to  the  problem  of  the  previous  an- 
nual report.  In  Figure  2-27  a typical  set  of  results  is  shown 
for  performance  as  a function  of  waviness  amplitude . The  param- 
eter values  used  are  the  same  as  those  in  Reference  [1]  for  base 
case  2.  The  combined  roughness  a is  taken  as  1/3  c where  c was 
one  half  of  the  total  roughness  height  in  the  previous  report. 

The  truncation  level  was  taken  at  1.5  ct^  or  tr  = 1.06.  The  re- 
sults for  the  truncated  roughness  model  can  be  compared  to  the 
original  base  case  2 results  reproduced  as  Figure  2-28.  For  the 
truncated  model,  the  percent  hydrodynamic  load  support  increases 
more  rapidly  and  reaches  100%.  This  behavior  is  expected  because 
as  stated  earlier,  the  use  of  the  Christensen  [20]  and  Tonder  [21] 
assumption  is  conservative,  i.e.,  greater  hydrodynamic  load  sup- 
port is  actually  achieved  than  predicted  by  this  model.  Leakage 
and  mean  film  thickness  are  plotted  as  divided  by  three  in  Figure 
2-27  so  as  to  be  directly  comparable  between  the  two  figures. 
These  results  are  similar  for  the  two  models.  Friction  coeffi- 
cient for  the  truncated  model  is  significantly  lower  and  was 
plotted  on  a different  scale  than  in  Figure  2-28  so  that  this  ef- 
fect could  be  shown. 

At  this  point  it  would  appear  that  the  results  based  on 
the  two  different  models  are  similar  in  most  ways  but  that  the 
truncated  model  provides  a larger  load  support.  One  question 
that  remains  is  the  choice  of  truncation  level.  Table  2-7  shows 


Table  2-7.  Truncation  Level  Study 


p = 5.55  x 10~4  p =■  5.55  x lo'2  p.  = p - 0.0 
o_  m_  c_ 


p * 0.1  p n =*  3 r.  =0.9  r /o  * 3.0  x 10^ 
s in  x o 

a a 


B = 0.75  * 02  - CT  h3(J  = 1.2 


t = 1.25  a. 

t = 1.50  a. 

V =*  1.75  a 

r 1 

r 1 

r 1 

% 

99.8 

89.5 

-P 

u 

o 

S'  V. 

u 

0.0019 

0.0124 

Q,  <0 

5 3 

2a 

0.1022 

0.1066 

4J  C 

g 2 

h 

o 

3.0555 

3.140 

o 

a 

the  effect  of  various  truncation  levels  for  a fixed  waviness.  The 
results  show  that  more  severe  truncation  favors  higher  hydrody- 
namic load  support.  At  less  severe  truncation  levels,  the  results 
show  that  the  seal  ring  cannot  be  supported.  These  results  suggest 
that  in  reality  the  truncation  level  must  adjust  in  some  manner  to 
the  loading  conditions. 

The  previous  results  suggest  that  useful  seal  modeling  may 
be  possible  using  even  the  simplest  of  interpretations  for  surface 
roughness.  However,  refined  models  will  require  a full  under- 
standing of  the  intricacies  discussed  herein.  What  is  needed  at 
this  point  is  a much  better  understanding  of  the  actual  surfaces 
which  are  developed  during  the  wear  process.  The  test  apparatus 
described  in  Chapter  4 will  be  used  to  further  this  understanding. 
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When  the  nature  of  worn  surfaces  is  better  understood  experi- 
mentally, then  a better  interpretation  of  the  preceding  results 
can  be  made. 


CHAPTER  3 


WAVY  SEAL  DESIGN 

PARAMETER  STUDY  RESULTS 

There  are  a number  of  design  parameters  that  have  a direct 
effect  on  the  performance  of  a wavy  mechanical  face  seal.  Among 
these  parameters  are  seed,  face  width,  surface  roughness,  number  of 
waves  around  the  circumference  of  the  seal,  and  the  amplitude  and 
shape  of  the  waves.  The  values  of  these  parameters  may  be  con- 
trolled to  some  degree  by  a face  seal  designer. 

As  a first  step  in  developing  a rational  design  procedure, 
the  effect  of  these  parameters  on  seal  performance  must  be  studied. 

The  following  studies  are  based  on  base  case  2 of  Reference  [1]  shown 
in  Table  3-1.  Although  these  results  were  presented  in  the  previous 
annual  report  [1],  they  will  be  reviewed  here  to  facilitate  discussion 
of  the  design  procedure. 

Pressure 

The  effect  of  changing  the  outside  pressure,  p , is  shown 

o 

in  Figure  3-1.  As  the  outside  pressure  Pq  is  increased,  the  load 
that  must  be  balanced  also  increases.  Consequently,  more  asperity 
contact  is  needed  and  hQ  decreases  accordingly  to  balance  the  load. 

The  fraction  of  the  load  supported  by  fluid  film  pressure  decreases 
accordingly.  The  leakage  increases  with  increasing  pressure  be- 
cause the  radial  pressure  gradient  is  increasing. 

The  coefficient  of  friction  remains  constant.  As  p is 

o 

increased,  the  total  fluid  friction  force  remains  essentially 


Table  3-1.  Base 

Case 

Parameters 

FACE  SEAL  APPLICATION  PARAMETERS 

Fluid  Viscosity 

n 

3.05xl0"4  Pa*s  (93°C  HO) 
(6. 37x10"®  lb* s/ft2)  2 

Outside  Pressure 

Po 

2.87  MPa  (416  psi) 

Angular  Speed  of  Shaft 

U) 

188  r/s  (1800  rpm) 

Mean  Radius 

r 

c 

47. 54  mm  (1.87  in. ) 

FACE  SEAL  DESIGN  PARAMETERS 

Fixed  Parameters 

Ring  Stiffness  Ratio 

A 

6.0 

Compressive  Strength 
of  Carbon 

Pm 

250  MPa  (36000  psi) 

Balance  Ratio 

B 

0.75 

Variable  Parameters 
(Base  Case  Shown) 

Surface  Roughness 

C 

0.5  um  (20  uin.) 

Face  Width 

Inside  Radius 

r . 

1 

• 45  mm  (1.77  in. ) 

Outside  Radius 

r 

o 

50  mm  (1.97  in.) 

Number  of  Waves 

n 

3 

Net  Waviness  Amplitude 

h3 

0.1  ^m  (4  uin.) 

Moment  of  Inertia 

I 

x 

4830  mm4  (0.0116  in.4) 

Modulus  of  Elasticity 

E 

21  GPa  (3.0x10s  psi) 

Non-Dimensional  Parameters 

Outside  Pressure 

P0 

0.005 

Compressive  Strength 

0.50 

Inside  Radius 

r . 

i 

0.9 

Net  Waviness  Amplitude 

S3 

0.2 

Ring  Stiffness 

K 

2. 5x10' 7 

Shear  Strength 
of  Asperities 

K 

0.05  (0.1  p ) 
m 

Friction  Parameter 

r /c 
o 

1.0x10s 

- 
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constant  and  the  coefficient  of  fluid  friction  decreases.  However, 

the  coefficient  of  mechanical  friction  increases  as  p is  increased. 

o 

The  addition  of  these  two  coefficients  of  friction  results  in  an  es- 
sentially constant  total  coefficient  of  friction. 

Face  Width 

The  results  of  changing  the  seal  face  width  are  shown  in 
Figure  3-2.  The  percent  hydrodynamic  load  support  increases  with 
increasing  face  width  as  more  area  becomes  available  for  load  sup- 
port. The  coefficient  of  friction  decreases  as  hydrodynamic  action 
improves.  The  mean  film  thickness  increases  as  less  penetration  is 
needed  to  develop  the  required  asperity  load  support.  The  leakage 
decreases  because  the  radial  pressure  gradient  is  decreasing  with 
increasing  face  width.  It  appears  that  a seal  with  a wide  face  has 
several  advantages. 

Surface  Roughness 

Figure  3-3  shows  the  effect  of  surface  roughness  on  face  seal 
performance.  Due  to  the  fact  that  the  non-dimens ionalization  of 
many  of  the  parameters  involves  the  parameter  "c,"  the  dimensional 
values  of  the  base  case  parameters  are  shown  on  the  figure.  Cavi- 
tation occurs  only  when  the  surface  roughness  is  less  than  1.0  ym 
(40  yin.).  The  percent  hydrodynamic  load  support  rises  sharply 
when  the  surface  roughness  becomes  less  than  0.64  urn  (25  yin.). 

One  hundred  percent  hydrodynamic  load  support  would  be  expected 
when  c * 0.33  ym  (13  yin.).  The  coefficient  of  friction  decreases 


considerably.  Due  to  the  non-dimensionalization  of  leakage,  the 


actual  leakage  decreases  more  than  the  curve  for  non-dimensional 


leakage  shows.  The  non-dimensional  mean  film  thickness  increases 
snarply.  However,  the  dimensional  mean  film  thickness  decreases 
as  the  surface  roughness  decreases.  With  regard  to  the  fraction 
of  hydrodynamic  load  support,  the  larger  roughness  offers  a much 
lower  resistance  to  fluid  flow  in  the  tangential  direction.  Hence, 
the  hydrodynamic  pressure  will  not  develop  to  as  great  an  extent 
as  for  the  lower  roughness  values.  The  coefficient  of  friction 
decreases  as  the  hydrodynamic  mechanism  improves  (with  lower  sur- 
face roughness  values) . The  leakage  will  increase  with  increasing 
surface  roughness  because  the  amount  of  area  available  for  leakage 
also  increases. 

Number~of  Waves 

The  effect  of  the  number  of  waves  around  the  face  seal  cir- 
cumference on  the  performance  parameters  is  illustrated  in  Figure 
3-4.  As  the  number  of  waves  increases,  the  percent  hydrodynamic 
load  support  increases.  This  is  due  to  the  fact  that  as  the  number 
of  waves  increases,  the  circumferential  wavelength  decreases.  With 
r ^ * constant,  the  bearing  width-to- length  ratio  increases  and  side 
leakage  is  reduced.  A greater  fraction  of  the  fluid  must  then  pass 
through  the  converging  film  region  and  hydrodynamic  pressure  becomes 
higher. 


As  the  hydrodynamic  forces  support  more  of  the  load,  less 
asperity  contact  is  needed  and  hQ  increases  accordingly.  The  leak- 
age remains  essentially  constant.  There  are  more  regions  through 


which  leakage  can  flow  (with  increasing  the  number  of  waves) , but 
each  individual  region  is  smaller. 

Wave  Amplitude 

The  effects  of  changing  the  net  waviness  amplitude  are  shown 
in  Figure  3-5.  The  hydrodynamic  load  support  increases  dramatically 
between  0.0  <_  h^  Z.  2.0.  However,  as  the  net  waviness  amplitude  con- 
tinues to  increase,  the  percent  hydrodynamic  load  support  curve  be- 
gins to  flatten  out.  This  suggests  that  there  is  some  limiting 
value  for  net  waviness  amplitude  beyond  which  the  hydrodynamic  ac- 
tion will  not  improve.  The  coefficient  of  friction  decreases  as 
the  hydrodynamic  lubrication  increases. 

The  mean  film  thickness  increases  as  the  net  waviness  ampli- 
tude increases.  With  the  hydrodynamic  mechanism  playing  a greater 
role  in  supporting  the  load,  less  asperity  contact  is  necessary. 

As  a result,  h increases.  As  h and  h„  increase,  the  leakage 
o 0 3 

also  rises  sharply. 

Wave  Shape 

All  of  the  previous  studies  herein  were  made  using  a simple 
cosine  wave  shape.  As  reported  previously  [1],  a study  was  made 
to  determine  if  other  wave  shapes  could  be  found  which  provided 
additional  load  support  while  maintaining  a constant  and  given 
leakage.  According  to  the  results  of  the  study,  the  cosine  wave 
shape  was  found  to  be  the  best  choice . No  other  wave  shape  tried 


improved  the  results. 


APPLICATION  OF  PARAMETER  STUDY  RESULTS  TO  FACE  SEAL  DESIGN 

So  far,  no  consideration  has  been  given  to  practical  values 
of  the  preceding  parameters  for  optimum  face  seal  design.  Each  of 
the  face  seal  parameters  is  now  considered  from  this  viewpoint. 

Face  Width 

The  results  of  the  face  width  study  indicate  that  better 
face  seal  performance  is  achieved  as  the  seal  face  width  is  in- 
creased. However,  the  heat  generated  at  the  seal  faces  is  approxi- 
mately proportional  to  the  seal  face  area.  As  the  heat  generation 
increases,  thermal  distortion  such  as  thermal  coning  increases. 
Therefore,  a larger  radial  taper  will  result  as  the  seal  face  width 
is  increased.  Radial  taper  effects  cause  uneven  seal  face  wear  (in 
the  radial  direction)  and  increased  leakage.  In  addition,  the  re- 
sulting increased  temperature  at  the  seal  interface  may  cause  vapor- 
ization of  the  liquid  between  the  seal  faces,  which  in  turn  may 
cause  the  condition  of  lubrication  between  the  seal  faces  to  change 
entirely.  As  a result,  the  present  practice  in  face  seal  design  of 
using  a seal  face  width  of  approximately  ri  0.9* appears  to  be  a 
logical  first  choice  for  a wavy  face  seal  design. 


Surface  Roughness 

Figure  3-3  illustrates  the  importance  of  surface  roughness 
on  the  performance  of  a mechanical  face  seal.  Considerable  hydro- 
dynamic  lubrication  can  be  achieved  if  the  asperity  roughness  height 
is  low.  For  the  numerical  study,  the  value  of  the  surface  roughness 
can  be  lowered  in  an  arbitrary  manner.  The  range  of  values  of  surface 

*For  seals  up  to  a few  inches  in  diameter. 
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roughness  which  can  be  obtained  in  practice  must  be  determined. 

Results  presented  previously  [1]  indicate  that  some  face 
seals  presently  have  surface  roughness  amplitudes  low  enough  to 
allow  hydrodynamic  lubrication.  However  Figure  3-3  shows  that 
even  lower  values  of  roughness  amplitude  are  desirable . Since 
little  research  has  been  conducted  on  the  relationship  between  ma- 
terial properties  and  surface  roughness,  it  must  be  assumed  at 
this  point  that  material  pairs  can  be  selected  that  result  in  a 
sufficiently  low  surface  roughness  so  that  a wavy  seal  design  can 
be  made  practical.  Results  from  the  experiments  discussed  in 
Chapter  4 will  provide  further  information  on  this  important  sub- 
ject. For  now  it  is  assumed  that  c = 20  uin.  is  attainable  in 
practice . 

Waviness  Amplitude 

Initial  waviness  is  that  amount  of  waviness  that  is  some- 
how imposed  upon  a mechanical  face  seal  ring  when  it  is  completely 
unloaded.  When  a face  seal  is  placed  into  operation,  the  sealed 
pressure  (closing  force)  flattens  the  wavy  face  seal  ring  against 
the  face  seal  seat.  As  a result,  the  initial  waviness  is  reduced 
by  some  amount  resulting  in  a net  waviness,  which  is  the  parameter 
used  in  previous  discussions.  Net  waviness  is  the  important  param- 
eter when  considering  the  waviness  effects  upon  hydrodynamic  lubri- 
cation. The  initial  waviness  is  important  also  because  it  serves 
as  the  source  of  waviness.  The  net  waviness  amplitude  is  related 
to  the  initial  waviness  amplitude  primarily  through  the  stiffness 
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factor  of  the  face  seal  ring. 

Relationships  between  initial  and  net 
previously  [1J . The  most  important  of  these 


waviness  were  developed 
relationships  are : 
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2m4  + n2 
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K - stiffness  constant  for  ring 
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rc  ra<iius  to  seal  ring  centroid 


A - ring  constant 
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n - number  of  waves  around  the  seal  face 
X'  5*l  • F°utier  of  the  deflected  vevie.es 

PV  Pfal  " POUrier  c°e«icient  of  the  pressure 

For  any  particular  face  seal  n A anH  v . 

n,  a,  and  K are  known  constants  and 

% “d  \ CSn  be  ^termined  by  use  of  Equation  (2-62)  in 
Reference  (1], 

The  net  v.vieess  sepUtude.  h,.  is  Mth-»ticaUy  defined 


as 


h.  + h 
la,  da 
1 1 


(3-3) 
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Using  the  above  relationships,  the  results  in  Figure  3-6 
were  obtained.  Using  Figure  3-6  in  conjunction  with  Figure  3-5 
gives  the  initial  waviness  amplitude  required  for  a desired  net 
waviness  amplitude.  The  larger  the  value  of  ring  stiffness  ic, 
the  less  initial  waviness  amplitude  is  required  to  achieve  a par- 
ticular net  waviness,  as  would  be  expected.  It  is  also  apparent 
that  for  a given  initial  waviness,  there  may  be  as  many  as  three 
different  operating  conditions  (net  waviness  values) . (See  points 
A,  B,  and  C for  K ■ 5.0  x 10  on  Figure  3-6).  The  determining 
factor  as  to  whether  there  are  three  possible  operating  conditions 
or  only  one  is  the  stiffness  of  the  ring.  The  less  stiff  rings 
have  three  possible  operating  conditions,  though  only  two  may  be 
shown  in  Figure  3-6.  Actually,  for  the  lower  values  of  ic  the  curves 
will  continue  to  rise  when  h^  > 5.0.  Therefore,  there  will  be  three 
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points  of  intersection.  The  stiffer  rings  of  Figure  3-6  have  only 
one  operating  condition  for  a particular  initial  waviness  as  shown. 

Using  methods  developed  and  presented  in  the  previous  annual 
report  [1],  it  was  determined  that  operating  points  in  Figure  3-6 
for  which  the  h^  versus  h3  curve  has  a positive  slope  are  stable. 
Where  the  slope  is  negative,  operation  is  unstable.  Thus  operation 
can  occur  at  points  A and  C,  but  not  at  point  B in  Figure  3-6. 

Figure  3-6  also  shows  that  the  amount  of  initial  waviness  must  be 
carefully  controlled  in  some  cases.  For  example , if  operation  is 
occurring  at  point  A and  a small  increase  in  initial  waviness  oc- 
curs, the  stable  operation  point  will  move  beyond  point  C and  high 
leakage  will  result. 

To  provide  additional  design  information  on  the  effect  of 
initial  waviness,  additional  studies  were  made  for  larger  numbers 
of  waves.  All  parameters  were  held  at  the  base  case  values  except 
the  number  of  waves  was  changed  from  three  to  six  and  to  nine.  The 
performance  curves  and  the  required  initial  waviness  curves  for 
n = 3 are  illustrated  in  Figures  3-5  and  3-6,  respectively.  These 
are  to  be  compared  to  Figures  3-7  and  3-8  for  n = 6,  and  Figures 
3-9  and  3-10  for  n = 9.  Comparing  the  three  sets  of  curves,  a few 
generalizations  cam  be  made.  The  performance  curves  of  Figures  3-5, 
3-7,  and  3-9  indicate  that  the  percent  hydrodynamic  load  support 
and  face  seal  leadcage  increase  to  a higher  level  at  a faster  rate 
with  increasing  net  waviness  amplitude  for  the  higher  values  of  n. 
Correspondingly,  the  coefficient  of  friction  decreases,  for  a par- 
ticular net  waviness,  as  n increases. 


Figure  3-9 . Effect  of  Net  Waviness  (n 


The  /required  initial  waviness  curves  for  n • 3,  6,  and  9 are 
shown  in  Figures  3-6,  3-8,  and  3-10,  respectively.  At  the  higher 
number  of  waves  the  multiple  operating  conditions  for  a particular 
initial  waviness  disappear.  That  is,  the  stability  problems  ap- 
parent for  n ■ 3,  Figure  3-6,  vanish  at  the  higher  values  of  n.  It 
appears  that  ary  stability  problem  can  be  avoided  by  imposing  a suf- 
ficient number  of  waves  on  the  seal  face. 

It  is  obvious  from  Figure  3-10  that  the  face  seal  ring  stiff- 
ness constant  K has  less  influence  on  the  initial  waviness  versus 
net  waviness  behavior  as  the  number  of  waves  is  increased.  The 
reason  for  this  is  that  the  effective  stiffness  of  the  seal  ring  in- 
creases greatly  as  n increases  (see  Equations  3-1  and  3-2)).  Con- 
sequently, the  face  seal  deflects  less  and  the  net  waviness  ampli- 
tude approaches  the  value  of  the  initial  waviness  as  a limit. 

Number  of  Waves 

The  results  of  the  parameter  study  concerning  the  number  of 
waves  on  the  seal  face  indicate  that  by  increasing  the  number  of 
waves,  face  seal  performance  improves.  Figures  3-7  and  3-9  show 
that  stability  problems  are  eliminated  by  using  a higher  number  of 
waves.  Whether  or  not  there  exists  a physical  limit  on  how  many 
waves  can  be  elastically  imposed  on  the  seal  face  must  now  be  de- 
termined. 

The  method  of  imposing  waviness  on  a seal  face  to  be  con- 
sidered here  is  the  elastic  method  discussed  in  Reference  [1], 
where  two  0-rings  provide  a bending  moment  when  pressure  is  applied. 
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The  reason  for  considering  an  elastic  method,  as  opposed  to  lapping 
the  waviness  into  the  face , is  that  with  such  a method  the  waviness 
can  be  moved  around  the  seal  ring,  thereby  evenly  distributing  the 
wear  on  the  seal  face  and  preserving  the  waviness  (see  Reference 

[ID  • 

Figure  3-11  shows  the  face  seal  ring  cross-section  and  the 
waviness  cylinder  (with  only  one  wave  shown) . The  wavy  lines  in 
Figure  3-11  are  the  outside  edges  of  the  O-ring  grooves.  The  area 
between  the  wavy  lines  is  to  be  pressurized  to  cause  an  angularly 
varying  bending  moment  about  the  seal  ring  centroid,  thereby  de- 
forming the  seal  ring  into  the  desired  waviness. 

The  amplitude  of  the  resulting  waviness  is  a function  of 
many  parameters  including  the  pressure  used  between  the  O-rings  and 
the  number  of  waves  around  the  face  seal.  Assuming  that  the  O-ring 
grooves  are  in  the  shape  of  a sine  function,  the  desired  wave  shape, 
the  moment  that  results  when  the  area  between  the  O-rings  is  pres- 
surized takes  the  form 

M = M sin  n0  (3-6) 

o 

The  maximum  moment  that  can  be  achieved  for  a given  pressure  is  of 
interest  since  any  lesser  moment  will  also  be  possible.  That  is, 
let  sin  n9  = 1.0.  From  Figure  3-11,  the  maximum  bending  moment  is 

- «)pb  (3-7) 

where  a - amplitude  of  waviness  of  the  O-ring, 
i - seal  ring  length, 
pfa  - pressure  used  to  cause  bending, 


(2a) 


<L 
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M - maximum  distributed  bending  moment  (in  lb.  per  in.  of 
o 

circumference) • 

It  would  appear  from  Equation  (3-7)  that  an  optimum  value 
of  "a"  exists  to  maximize  Mq.  Since  t and  pfa  are  constants,  the 

optimum  value  of  "a"  is  arrived  at 
I dM 

— 2.  = Pw  (1  - 4a)  =0  (3- 

da  b 

Therefore,  a = 1/ 4 is  the  optimum  value  of  "a."  Substituting  this 
value  for  "a"  back,  into  Equation  (3-7) 

m = ~ £2Pb  (3‘ 

o 8 D 

max 

Lebeck  [22]  has  derived  an  expression  that  relates  the 
initial  waviness  amplitude  h^  to  the  bending  moment  Mq. 


h.EI  (n2  - l)2 
i x 

M = ; 

° r3  (1  + A) 

c 


(3-10) 


By  equating  the  two  expressions  for  Mo  (Equations  3-9)  and 
(3-10)),  the  desired  relationship  among  hi#  Pb,  and  n is  obtained. 


£2r3  ( 1 + A) 
c 

\ “ 2 : 

1 8EI  (n*4  - 1) 
x 


(3-11) 


Non-dimensionalizing  Equation  (3-11),  the  following  expression  is 


obtained: 


o 


T pb 


(3-12) 


For  a mechanical  face  seal  to  operate  in  a particular,  de- 
sired manner  (region),  there  exists  a value  of  hn,  net  waviness. 
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that  must  be  present.  To  obtain  this  h , an  initial  waviness  h. 

n i 

must  be  imposed  on  the  seal  face,  defined  by  Equation  (3-12).  If 

the  waviness  pressure  5,  is  held  constant,  then  for  any  specific 

o 

value  of  n (all  other  parameters  held  constant)  there  is  a maximum 

that  can  be  imposed.  Any  lesser  value  of  can  bfe  achieved. 

That  is,  there  exists  a range  of  h. , and  thus  h , that  is  possible, 

a n 

given  pfa,  K,  n,  and  the  other  face  seal  parameters  in  expression 
(3-12).  For  example,  for 


t_ 

r 

o 


0.5 


A * 6.0 


p * p a 0.005 
b o 


n - 3,  6,  9 


(3-13) 


K = variable 

the  following  shown  in  Table  3-2  are  possible. 

A comparison  is  made  between  the  maximum  h^  possible  for  a 

specific  case  and  the  iv  required  to  result  in  a particular,  desired 

h . Referring  to  Figure  3-6  , all  stable  values  of  h , within  the 

n J 

range  of  the  figure,  are  possible  for  each  value  of  K. 

However,  a comparison  between  the  values  in  Table  3-2  for 
n =*  6 and  n * 9 and  Figures  3-8  and  3-10  shows  that  not  all  opera- 
ting conditions  are  possible  (for  the  example  face  seal) . For 

n * 6 and  K » 1.25  * 10  ^ , all  desired  h.  shown  in  Figure  3-8  can 

o 

be  achieved.  However,  for  all  other  K curves  of  n * 6,  and  for 
all  K curves  for  n * 9,  there  are  intersections  of  the  h 
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and  the  h.  versus  h curves.  The  physically  possible  operating 
i n 

conditions  are  those  to  the  left  of  the  point  of  intersection  for 

-7 

each  curve.  For  example,  with  n = 6 and  K = 2.5  x 10  there  is 

a maximum  possible  initial  waviness  that  can  be  achieved  of  3.6 

non-dimensional  units.  Correspondingly,  referring  to  Figure  3-8, 

the  maximum  h,  that  can  result  is  h.  = 2.7  units.  Therefore, 

6 o 

0 < h.  < 2.7  is  the  range  of  possible  operating  conditions  for 
— 6 

this  particular  case.  Again  referring  to  Figure  3-8,  as  the  seal 

ring  stiffness  increases,  the  amount  of  initial  wavir.ess  that  can 

be  imposed  elastically  decreases,  as  does  the  resulting  h&  that 

can  be  achieved.  In  fact,  h,.  reaches  such  a low  value  that  it  is 

6 

doubtful  that  any  appreciable  hydrodynamic  action  will  occur. 

The  same  behavior  exists  for  n = 9 (Figure  3-10) . The  ranges 
of  possible  operating  conditions  are  to  the  left  of  the  points  of 
intersection.  Again,  the  range  of  possible  hg  decreases  rapidly 
with  increasing  ic,  to  such  a point  that  hydrodynamic  action  becomes 


very  small. 


Table  3-3  summarizes  the  maximum  net  waviness  amplitudes 


that  result  when  the  maximum  initial  waviness  amplitudes  for  stable 

operating  conditions  are  imposed.  This  information  is  obtained 

from  Figures  3-6,  3-8,  and  3-10. 

Table  3-3.  Maximum  Met  Waviness  Amplitudes 
for  n * 3,  6,  9 
(for  stable  operation) 


Ring  Stiffness  (K) 

hw 

^max 

3 max 

1.25  x 10~7 

0.75 

>5.00 

1.00 

2.50  x 10~7 

0.75 

2.70 

0.42 

5.00  x 10~7 

0.75 

0.18 

1.00  x 10“ 5 

> 5.00 

mm 

0.08 

2.00  x 10_s 

> 5.00 

0.10 

0.05 

_ 

It  is  important  to  note  that  the  behavior  of  h^  (shown  in 

max 

Table  3-3)  results  from  the  instability  problems  associated  with 
the  n = 3 case,  displayed  graphically  in  Figure  3-6,  and  discussed 
in  detail  previously. 

It  can  be  seen  from  the  performance  curves  (Figures  3-5, 

3-7,  and  3-9)  that  as  the  net  waviness  amplitude  increases,  the 

face  seal  leakage  and  hydrodynamic  load  support  also  increase. 

Therefore,  since  the  h maximum  is  known  for  each  value  of  K and 

n 

n studied  (see  Table  3-3),  the  maximum  face  seal  leakage  and  maxi- 
mum hydrodynamic  load  support  are  also  known.  This  information  is 
displayed  graphically  in  Figures  3-12  and  3-13.  The  curves  of 


Figure  3-12  connect  points  of  maximum  percent  hydrodynamic  load 


is 
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support  for  n = 3,  6,  and  9,  and  1.25  *10  K 0.0  x 10 
Figure  3-13  is  the  same  type  of  graph  as  Figure  3-12  only  for  maxi- 
mum face  seal  leakage.  The  possible  operating  ranges  are  to  the 
left  of  the  n = 6 and  n = 9 curves  and  below  the  n = 3 curves. 

The  step-like  behavior  of  the  n = 3 curves  graphically  por- 
trays the  instability  problems  associated  with  this  case  (see 
Table  3-3)  . In  addition,  Figures  3-12  and  3-13  emphasize  the  fact 
that  as  n is  increased,  the  range  of  possible  operating  conditions 
is  reduced.  With  these  figures  one  can  tell  at  a glance  the  physi- 
cally possible  hydrodynamic  load  support  and  leakage  ranges  as  a 
function  of  K and  n. 

EXAMPLE 

In  order  to  illustrate  how  the  previous  information  can  be 
used  in  the  design  of  a wavy  mechanical  face  seal,  an  example  is 
presented.  Two  design  procedures  are  used  for  comparison  and  il- 
lustration purposes. 

The  application  that  is  considered  for  this  example  is  that 
of  the  "base  case"  (see  Table  3-1) . The  maximum  allowable  leakage  for 
this  example  is  taken  as  2 cc/min.  This  value  for  leakage  is  chosen 
because  it  offers  the  greatest  number  of  possible  design  complica- 
tions including  stability  considerations,  and  because  this  value 
represents  an  acceptable  leakage  for  many  applications.  Non- 
dimensional  leakage,  Q,  is  0.0878. 
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There  are  a few  assumptions  that  must  be  made  concerning  the 
mechanical  face  seal  application,  materials,  and  face  seal  geometry. 


1)  The  seal  face  width  is  r.  1 0.9. 

2)  The  surface  roughness  of  c » 0.5  um  (20  uin.)  is 
attainable . 

3)  The  traveling  wave  concept  is  to  be  used  in  order  that 
the  waviness  will  not  wear  away. 

Three  of  the  five  design  parameters  studied  in  this  investi- 
gation have  been  fixed  at  certain  values  for  this  example.  The 
reasons  for  this  were  di  cussed  previously.  The  two  design  param- 
eters that  are  still  to  be  examined  are  the  number  of  waves  around 
the  seal  face,  and  the  amplitude  of  the  net  waviness.  The  objective 
is  to  attain  the  maximum  amount  of  hydrodynamic  load  support  while 
maintaining  the  leakage  at  or  below  the  maximum  allowable  leakage. 

Procedure  1 

The  most  useful  figures  are  3-5  through  3-10.  These  figures 
relate  the  two  remaining  design  parameters.  From  Figure  3-5  for 
n * 3,  the  h3  that  gives  Q <_  0.0878  is  <_  O-SO.1  From  Figure 
3-6  for  h3  * 0.60  and  K * 2.5  * 10  7 (the  base  case  value  of  K) , 
the  stable  h^  required  is  h^  = 36.8  units.  For  n = 3 all  values 
of  initial  waviness  for  Figure  3-6  are  possible.  Although  an 


The  figures  actually  require  log  Q.  Therefore,  log  Qmax  = -1.075. 
Figures_3-5,  3-7,  and  3-9  indicate  that  the  largest  h^  that  still 
limits  Qmax  <_  0.0878  and  provides  the  highest  percent  hydrodynamic 
load  support  is  best,  so  the  inequality  will  no  longer  be  considered. 


= 36.8  units  is  physically  possible  to  achieve,  this  particular 
value  of  h^  is  very  close  to  the  peak  of  the  h.  versus  h^  curve  for 
K = 2.5  * 10  7 (see  Figure  3-6).  Consequently,  if  cannot  be 
imposed  precisely  enough,  stability  problems  are  likely  to  occur. 
This  would  not  be  a recommended  operating  condition.  However,  for 
the  purpose  of  illustration,  it  is  assumed  that  h^  can  be  made  36.8 
units.  Therefore,  a design  solution  for  this  example  is  n = 3, 

= 0.60  and  h^  = 36.8,  with  all  other  design  parameters  as  pre- 
viously specified.  Referring  to  Figure  3-5,  the  percent  hydro- 
dynamic  load  support  that  is  expected  for  this  design  is  84.5  per- 
cent, a significant  gain  over  the  lower  limiting  value  of  69  per- 
cent. The  coefficient  of  friction  is  0.0178.  This  result  is 
summarized  in  Table  3-4. 

Table  3-4.  Example  Face  Seal  Design  Results 

(Procedure  1) (Q  = 0.0878) 
max 


n = 3 

n = 6 

n = 9 

n = 9 (modified) 

h. 

36.3 

1.98 

0.68 

0.775 

h 

n 

0.60 

0.55 

0.42 

0.525 

% 

84.5 

95.3 

94.50 

98.20 

u 

0.0178 

0.0065 

0.0066 

0.CO36 

The  same  procedure  used  for  n = 3 is  used  for  n * 6.  The 
results  are  shown  in  Table  3-4.  The  h^  required  is  physically  pos- 
sible. There  is  a significant  increase  in  percent  hydrodynamic 
load  support  and  a corresponding  decrease  in  the  coefficient  of 
friction.  There  is  no  stability  problem  for  this  design. 
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Finally,  the  design  procedure  is  used  for  n = 9.  This  case 
is  different  from  the  n = 3 and  n = 6 cases  and  is  explained  in  de- 
tail. From  Figure  3-9  for  log  =*  -1.057,  the  hg  required  is 

hg  » 0.525.  Referring  to  Figure  3-10  for  hg  = 0.525,  the  h^  re- 
quired is  h.  * 0.75.  However,  this  value  is  larger  than  the  maxi- 
mum h^  that  can  be  achieved  (see  Table  3-2) . Consequently,  the 
highest  h^  that  cam  be  attained  (for  this  example)  is  h^  = 0.68. 

The  hg  that  corresponds  to  h^  = 0.68  is  hg  =»  0.42.  Referring  to 
Figure  3-9  for  hg  * 0.42,  the  percent  hydrodynamic  load  support 
is  % = 94.50  percent,  while  the  coefficient  of  friction  is 
U ■ 0.0066.  Table  3-4  summarizes  the  results  of  this  example 
face  seal  design. 

So  far  in  this  example  K has  been  restricted  to  the  base 
case  value . In  reality  n and  K would  be  allowed  to  vary  through 
a wider,  more  precise  range.  If  for  n • 9,  K were  allowed  to  as- 
sume any  value  and  in  particular  K = 2.33  * 10  7,  the  results  ob- 
tained make  the  n • 9 design  more  desirable.  (See  Table  3-4  for 
n = 9 (modified) ) . 


Procedure  2 

Another,  more  simplified  design  procedure  can  also  be  used 
for  this  face  seal  example  design.  Figures  3-12  and  3-13  contain 
all  of  the  most  important  information  on  Figures  3-5  through  3-10 
in  a very  concise  form.  In  addition,  K can  assume  any  value  in 
the  range  of  K that  has  been  examined. 

Referring  to  Figure  3-13,  the  points  of  intersection  of  the 
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leakage  curves  for  n = 3,  6,  and  9,  and  log  Q = -1.057  are  de- 

max 

termined.  (For  n = 3,  since  there  is  no  point  of  intersection,  an 
approximation  will  have  to  be  made  throughout  the  design  procedure.) 
These  points  of  intersection  define  the  optimum  ic  for  each  value 
of  n.  Knowing  the  optimum  values  of  K,  Figure  3-12  is  used  to 
determine  the  maximum  percent  hydrodynamic  load  support  that  can 
be  obtained  for  each  value  of  n.  The  results  are  shown  in  Table 
3-5. 

Table  3-5.  Example  Face  Seal  Design  Results 
(Procedure  2) (Qmax  = 0.0878) . 


■pn 

Ball 

n =*  6 

n * 9 

K (optimum) 

Determined  by 

7.55  x 10-7 

2.19  x lCf7 

Stability 

% 

<87.2 

86.7 

95.0 

The  results  shown  in  Table  3-5  differ  from  those  in  Table 
3-4  for  the  following  reasons.  Error  has  been  introduced  in  the 
graphical  method  used  in  producing  the  figures,  particularly  Fig- 
ures 3-12  and  3-13.  If  all  of  the  curves  shown  in  all  of  the 
figures  were  continuous,  the  results  in  Table  3-4  and  3-5  would 
be  identical.  However,  the  data  points  are  quite  widely  spread 
and  are  linearly  interpolated  causing  considerable  error. 

The  major  advantage  of  Procedure  2 over  Procedure  1 is  its 
speed  and  simplicity.  However,  an  important  disadvantage  of  Pro- 
cedure 2 is  demonstrated  with  the  n = 3 case.  For  face  seals  that 


131 


have  a stability  problem  (no  point  of  intersection  for  the  Q * 

max 

constant  curve  and  the  leakage  curve),  there  is  not  enough  informa- 
tion contained  in  Figures  3-12  and  3-13  to  determine  the  practical 
operating  conditions.  In  such  cases,  the  designer  must  use  the 
curves  of  Procedure  1.  In  addition,  although  Procedure  2 indicates 
the  possible  operating  conditions  for  a face  seal  design,  h.  must 
still  be  determined  from  the  figures  used  in  Procedure  1. 

In  summary,  the  wavy  face  seal  design  results  in  a consid- 
erable improvement  in  hydrodynamic  load  support  for  a given  leakage 
over  a non-wavy  face  seal.  Theoretically,  the  wear  rate  is  reduced 
by  a factor  of  17:1  for  the  n=9  case  where  % = 98.20  when  compared 
to  a flat  face  seal. 
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CHAPTER  4 


SEAL  TEST  APPARATUS 

To  date,  only  very  limited  studies  have  been  made  of  wavy 
seals  operating  in  water.  To  provide  additional  experimental  veri- 
fication of  the  model  described  herein  and  to  provide  facilities  for 
further  studies  related  to  other  seal  problems,  the  design  and  con- 
struction of  a mechanical  seal  test  apparatus  has  been  completed  as 
a part  of  this  project. 

BASIC  REQUIREMENTS 

A test  apparatus  designed  for  research  must  meet  more  re- 
quirements than  test  machines  used  for  study  of  commercial  seals. 
First  it  is  necessary  that  a wide  range  of  speed  and  sealed  pres- 
sure be  available  in  order  to  verify  modeled  behavior  over  a wide 
range  of  conditions.  Certain  otherwise  uncontrolled  variables 
such  as  shaft  runout,  seal  waviness,  and  seal  radial  taper  must  be 
very  carefully  controlled.  For  testing  the  lubrication  model  of 
interest,  it  is  necessary  that  a controllable  waviness  be  applied 
to  the  seal  faces,  and  in  order  to  reduce  the  effects  of  wear  this 
waviness  must  be  moved  circumferentially  with  respect  to  the  primary 
ring.  Leakage  from  the  seal  must  be  measured.  Torque  must  be 
measured  in  a manner  such  that  the  seal  torque  is  measured  inde- 
pendently of  bearing  torque  or  torque  arising  from  a second  seal 
used  in  a tandem  arrangement.  Many  of  the  distortions  commonly 
found  in  commercial  seals  must  be  minimized. 


To  meet  all  of  these  requirements,  the  design  of  such  a test 
apparatus  must  depart  somewhat  from  that  found  in  commercial  prac- 
tice. On  the  other  hand,  it  is  important  that  the  apparatus  be 
capable  of  simulating  commonly  found  seal  installations.  In  fact, 
it  is  useful  to  be  able  to  install  a commercially  available  stand- 
ard seal  in  the  test  rig. 


SPECIFICATIONS 


The  test  apparatus  specifications  which  satisfy  the  above 
basic  requirement  are  set  forth  below. 

Configuration:  Horizontal  shaft,  single  seal  with  torque  measured 
on  the  seal  housing.  Rotating  secondary  ring. 

Siz® : Accepts  up  to  3-5/8  inch  nominal  shaft  size  seal. 

4 inch  nominal  seal  face  diameter. 


Speed : 
Pressure : 
Bearings : 


Motor: 
Fluid: 
Balance : 
Seal  Rings: 

Waviness : 


500-4000  RPM.  9-70  ft/sec  for  a 4 inch  seal. 

0-1000  psi  continuous.  Vessel  is  rated  at  2000  psi. 
Five  ABEC-7  bearings  on  precision  spindle.  Bearings 
are  rated  at  a maximum  of  12000  lb  thrust  at  4000 
RPM  (corresponds  to  a 4 inch  seal  at  1000  psi) . 

5 Hp  variable  speed  drive  (belt  type) . 

Materials  are  selected  for  seawater  service . 

Balance  ratio  of  1.0  (lower  values  can  be  used). 

Zero  moment  design.  Designed  to  minimize  the  ef- 
fects of  pressure  caused  radial  taper. 

Seal  waviness  is  controllable  by  variable  ring 
distortion.  It  can  be  moved  relative  to  the 
primary  ring. 


Temperature : 


I 


Temperature:  Can  be  controlled  at  any  given  temperature  above 

100 °F. 

Quantities  Speed,  pressure,  temperature,  waviness,  torque. 

Measured: 

and  leakage. 


DESIGN 

Figure  4-1  shows  the  general  layout  of  the  test  apparatus. 
The  bolted  vessel  on  the  right  hand  end  contains  the  pressurized 
seal.  Immediately  to  the  left  is  the  waviness  drive  plate  and  the 
bird  cage-like  torque  transducer.  The  entire  test  vessel  is  sup- 
ported cantilever  style  by  the  torque  transducer.  Further  to  the 
left  is  the  bearing  housing  which  is  attached  to  the  frame.  The 
variable  speed  motor  drive  is  located  to  the  left  of  the  bearing 
housing . 

On  the  level  below  the  test  vessel  in  Figure  4-1  is  located 
all  of  the  auxiliary  equipment.  To  the  left  end  is  the  oil  supply 
system  including  a tank,  pump,  filter,  and  cooler.  On  the  right 
end  is  the  high  pressure  water  circulation  pump  heat  exchanger. 
Beyond  the  machine  to  the  right  is  the  control  panel  and  a rack 
containing  torque  and  temperature  recorders. 

Figure  4-2  shows  the  essential  features  of  the  seal  test 
apparatus  in  detail.  The  primary  and  secondary  seal  rings  are 


shown  as  parts  (27)  and 


Shaft  (l)  turns  on  five  precision 


angular  contact  bearings  designed  to  operate  under  the  large  thrust 
load.  Cooling  water  is  circulated  through  the  test  chamber  (^) 
through  openings  not  shown.  The  torque  tranducer  ^l)  acts  as 
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Fiyure  4-1 . Mechanical  Face  Seal  Test  Apparatus 


Figure  4-2.  Seal  Test  Apparatus  Assembly 


the  support  for  the  vessel.  The  pressure  vessel  is  designed  to  be 
readily  taken  apart  by  removal  of  bolts  (^3)  . The  secondary  ring 
support  rotor  is  removed  by  removing  capscrew  (^0)  . To  minimize 


axial  runout  at  the  seal  face , the  end  of  shaft  (T)  is  ground  true 
with  the  bearings.  Rotor  is  ground  both  where  it  mates  with 

the  shaft  and  secondary  seal  ring.  The  two  secondary  seal  ring 
faces  are  ground  parallel. 

Waviness  Drive 

The  traveling  wave  concept  is  shown  in  Figure  4-3.  The 
primary  ring  is  distorted  so  as  to  cause  waviness  on  the  face  of  the 
rings  when  mated  with  the  secondary  ring  or  seat.  After  a period  of 
time  this  wave  is  moved  relative  to  the  primary  ring  so  that  wear 
is  distributed  evenly  around  the  entire  ring  and  so  that  the  wave 
is  not  simply  worn  away.  Using  this  concept,  the  effects  of  wavi- 
ness in  contacting  face  seals  can  be  studied. 

Figure  4-4  shows  how  this  concept  was  implemented  in  the 
seal  test  apparatus.  The  drive  ring  is  designed  so  that  the  pri- 
mary ring  can  float  and  align  itself  but  cannot  rotate  about  the 
shaft  axis.  The  primary  ring  thus  takes  its  alignment  from  the 
face  of  the  rotating  secondary  ring.  The  waviness  drive  cylinder 
is  driven  by  a worm  at  low  speed  (2  revolutions  per  24  hours) . The 
waviness  drive  cylinder  seals  into  the  housing  with  two  O-rings  de- 
signed so  that  gas  pressure  can  be  supplied  through  the  end  plate 
and  into  the  waviness  cylinder.  Each  of  the  six  holes  on  the 
waviness  drive  cylinder  is  connected  to  this  gas  pressure . When 

• I 
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DISTORTED  RING  UNDER  PRESSURE 
AFTER  WAVE  HAS  MOVED 


Figure  4-3.  Traveling  Wave  Concept 
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Lfcfr'T  END  PlAl’E 


Figure  4-4.  Seal  Assembly 


the  primary  ring  is  properly  positioned  over  these  six  holes  and 
sealed  off  by  six  O-rings,  then  the  primary  ring  becomes  subjected 
to  six  small  regions  of  gas  pressure  on  its  inside.  Waviness  is 
caused  by  this  gas  pressure , because  of  the  staggered  location 
pattern  of  the  six  holes.  As  the  waviness  cylinder  is  turned  by 
the  worm,  the  waviness  pattern  moves  relative  to  the  seal  primary 
ring.  Figure  4-2  shows  additional  details  of  the  waviness  drive. 

Also  shown  in  Figure  4-4  is  the  secondary  O-ring  seal  be- 
tween the  waviness  cylinder  and  the  primary  ring.  This  O-ring 
seals  off  the  gap  between  these  two  members  while  allowing  the  pri- 
mary ring  to  float  to  align  itself. 

Primary  Ring 

The  primary  ring  is  shown  in  Figure  4-5.  The  material  se- 
lected for  initial  tests  is  Pure  Carbon  P658RC.  The  unusual  pro- 
portions of  the  ring  result  from  adjusting  the  geometry  so  that  no 
rotation  about  a circumferential  axis  is  caused  by  pressure  variation. 
The  primary  ring  is  driven  by  lugs  which  engage  two  notches.  Two 
as  opposed  to  a greater  number  of  notches  were  used  because  this 
arrangement  gives  the  primary  ring  maximum  freedom  to  float  radially. 
Second  harmonic  waviness  is  minimized  by  driving  the  ring  through 
the  centroid  of  its  cross  section.  Any  remaining  second  harmonic 
waviness  produced  by  drive  forces  will  be  easily  flattened  com- 
pared to  the  third  harmonic  waviness  forced  into  the  ring,  so  is  not 
considered  important.  The  secondary  seal  is  located  at  the  left  end 
on  the  inside  diameter.  The  balance  ratio  for  the  design  shown  is 
unity.  Springs  to  the  left  of  the  seal  provide  preload. 
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Secondary  Ring 


The  secondary  seal  ring  is  shown  in  Figure  4-6.  Rings  are 
fabricated  from  tungsten  carbide  and  ceramic.  This  ring  is  also 
of  a zero  pressure  moment  design.  An  unusual  feature  of  the  par- 
ticular design  is  that  the  mechanical  force  at  the  right  hand  side 
of  the  seal  has  been  reduced  to  nearly  zero  by  placing  the  O-ring 
seal  on  the  right  hand  face  as  opposed  to  the  outside  diameter 
(see  Figure  4-2) . It  is  useful  to  minimize  the  axial  force  to  re- 
duce coning  effects  caused  by  friction  at  the  load  bearing  point 
in  conjunction  with  radial  deformation  due  to  pressure  and  tempera- 
ture. The  secondary  seal  ring  is  held  in  place  by  spring  clips 
(not  shown)  which  engage  the  notches  shown  in  Figure  4-6. 

Torque  Transducer 

The  torque  transducer  is  shown  in  Figure  4-7.  The  diffi- 
culty in  designing  such  a torque  transducer  is  that  torque  to  be 
measured  is  small  whereas  the  axial  load  due  to  the  thrust  of  the 
seal  is  large.  Several  design  concepts  were  developed.  Figure 
4-7  shows  what  appears  to  be  the  best  compromise.  The  0.120  * 1 
inch  members  carry  the  axial  load  in  tension.  Due  to  their  small 
moment  of  inertia  about  a radial  axis,  these  members  have  in  com- 
bination a low  torsional  stiffness.  The  torsion  load  is  sensed 
by  the  two  reduced  section  beams.  Strain  gages  are  attached  to 
these  members  as  shown.  The  axial  load  effect  is  isolated  from 
these  members  by  the  cutout  shown.  The  torque  transducer  does 
have  some  small  amount  of  sensitivity  to  axial  loads  because  the 
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0.120  x 1 members  do  attempt  to  straighten  (when  twisted)  due  to 


the  axial  load.  This  effect  is  largely  offset  by  the  fact  that  the 
torque  sensing  beams  have  a very  high  torsional  stiffness.  Thus 
the  actual  torsional  deformation  of  the  unit  is  quite  small. 

The  circuit  diagram  for  the  torque  transducer  is  shown  in 
Figure  4-8.  A full  strain  gage  bridge  is  used  for  maximum  sen- 
sitivity and  best  thermal  compensation.  This  is  quite  important 
because  the  torque  transducer  is  mounted  to  the  bearing  housing 
which  reaches  a temperature  of  150  °F.  Test  results  show  up  to  a 7.5 
in.  lb.  zero  shift  caused  by  temperature  changes  of  the  torque  trans- 
ducer. Torque  readings  are  adjusted  for  this  shift. 

The  high  pressure  cooling  feed  to  the  vessel  (see  Figure  4-1) 
is  arranged  so  that  the  feed  tubes  sure  in  radial  planes  to  the 
machine  axis.  This  reduces  the  effect  of  pressure  on  the  zero  torque 
point  to  a minimum.  Tests  so  far  have  indicated  no  measurable  pres- 
sure effect  on  the  zero  point. 

Water  and  Gas  Systems 

Figure  4-9  shows  the  water  and  gas  flows  for  the  system. 
Entering  water  passes  through  valve  V2.  V2  is  electrically  con- 
trolled to  shut  off  in  the  event  of  major  leaks  in  the  system.  On 
the  low  pressure  water  side  of  the  system,  cooling  water  is  regu- 
lated by  servo  valve  MV1  and  then  passes  through  the  low  pressure 
side  of  heat  exchanger  HX1  and  then  returned  to  a drain. 

On  the  high  pressure  water  side,  water  is  pressurized  by 
air  driven  pump  PI  which  is  controlled  by  electrically  operated 
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Figure  4-8.  Torque  Transducer  Circuit 


Water  and  Gas  Systems 


Figure  4-9  (continued)  Water  and  Gas  Systems 


VI  - High  Pressure  Water  Pump  Electric  Solenoid  Air  Valve, 
Dayton  6X080,  1/4" 

V2  - Leak  Detector  Water  Shutoff  Solenoid  Valve,  Dayton 
6X082,  1/2" 

V3  - Water  Supply  Valve , 1/4"  Gate 

V4  - Water  Supply  Valve,  1/2"  Gate 

V5  - High  Pressure  Drain,  1/4",  3000  psi,  Nupro  Stainless 

V6  - High  Pressure  Bleed,  1/4",  3000  psi,  Nupro  Stainless 

V7  - Waviness  Gas  Supply,  1/4  Needle,  Alkron 
V8  - Accumulator  Valve 

V9  - Accumulator  Bleed  Valve,  1/4  Needle,  Alkron 

PI  - Air  Pump  - Hydraulic  Engineering  Co.,  SC10-500-16 

P2  - Circulation  Pump  w/High  Pressure  Inlet 

Centrifugal  type 

P3  - Sump  Pump  Centrifugal  Type  AMXCO  1/60  HP 

AC  - Accumulator,  Bladder  Type  Hydrill  K 1/4-3000  w/K31 
plastic  coating 

PG1  - 0-2000  psi  Gage 

PG2  - 0-2000  psi  Gage 

MV1  - Temperature  Control  Servo  Valve 
Honeywell  V5011A-1106  1/2" 

HX1  - Heat  Exchanger  Heliflow  4-X6C-10 
tube  and  shell  type 

R - Gas  Regulator  CGA580 

N - Nitrogen  Supply  Bottle 

SV1  - Blow  Out  Safety  Valve  Nupro  1200  psi  diaphragm 

S4  - High  Pressure  Controller,  100-1000  psi 
Honeywell  P634A 
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air  valve  VI  which  in  turn  is  controlled  by  pressure  controller  S4. 
Accumulator  AC  is  charged  with  nitrogen  to  75%  of  operating  pres- 
sure. This  reduces  pressure  surges  on  the  system.  This  high  pres- 
sure supply  system  is  required  to  only  make  up  water  that  has  been 
leaked  out  of  the  system  by  the  pump  seal  or  the  test  seal. 

The  high  pressure  water  supply  is  connected  to  the  cooling 
water  circulation  loop  as  shown.  Cooling  water  passes  through  the 
tube  side  of  HXl  and  is  circulated  through  the  test  vessel  by  cir- 
culation pump  P2.  Leakage  and  vessel  drains  are  collected  in  a 
sump  and  pumped  back  to  a drain. 

SV1  is  a burst  diaphragm  type  of  safety  valve  to  protect 
the  system  from  overpressure  in  the  event  of  a control  failure . 

High  pressure  nitrogen  gas  is  supplied  by  a bottle  and  reg- 
ulator to  either  charge  the  accumulator  or  to  pressurize  the 
waviness  cylinder. 

High  Pressure  Circulation  Pump 

The  high  pressure  circulation  pump  was  designed  and  built 
especially  for  the  project.  This  was  done  because  of  the  unavail- 
ability of  centrifugal  pumps  that  can  withstand  a high  inlet 
pressure . The  distinctive  feature  of  the  circulation  pump  shown 
in  Figure  4-10  is  its  heavy  walled  316  stainless  steel  housing. 

Oil  System 

Figure  4-11  shows  the  oil  lubrication  system  which  includes 
a pump,  filter,  air  cooled  heat  exchanger,  and  tank.  Switch  S2  shuts 
down  the  entire  test  apparatus  in  the  event  of  low  oil  pressure . 


-10.  High  Pressure  Circulation  Pump  Assembly 


-11.  Oil  Lubrication  System 


Control  Cabinet 


Figure  4-12  shows  the  circuitry  designed  for  the  remote  con- 
trol cabinet.  Switches  S6  through  S10  actuate  the  various  motors 
and  pumps  on  the  machine.  The  coil  on  REL3  is  connected  in  series 
with  several  important  safety  devices.  This  coil  is  deenergized 
in  the  event  of  oil  under  pressure  (Switch  S2  on  Figure  4-13) , 
oil  over  temperature  (CTL1) , seal  over  temperature  (CTL2) , or  a 
water  leak  (REL4).  When  REL3  is  deenergized,  all  systems  (all 
motors  and  pumps)  are  shut  off.  CTL3  is  the  controller  for  MV1 
in  Figure  4-9.  Seal  cooling  water,  seal  face,  and  oil  temperature 
are  monitored  by  thermocouples  connected  to  the  three  controllers. 
The  coil  of  REL4  is  connected  to  several  salt  soaked  and  dried 
probes  placed  around  the  machine.  The  probes  become  conductive 
when  contacted  with  water,  so  the  circuit  shuts  the  machine  down 
in  the  event  of  a major  water  leak. 

Machine  Electrical  Circuit 

Figure  4-13  shows  the  various  electrical  components  which 
are  mounted  on  the  machine  itself.  The  larger  motors  are  operated 
through  motor  starters  REL5,  REL6,  and  REL7,  which  include  over- 
load protection. 

SEAL  TEST  PROCEDURE 

The  following  procedure  for  using  the  apparatus  and  testing 
a seal  has  been  established. 

1.  Warm  up  torque  transducer  for  about  two  hours.  With 

secondary  seal  ring  removed,  the  bolt- tightening  bracket 
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Figure  4-12  Control  Cabinet  Circuit 


Figure  4-12  (continued)  Control  Cabinet  Components 
SI  - Main  Switch,  SPST  key  operated 

S3  - Low  Oil  Pressure  Override  Switch  - SPST  spring  toggle 

55  - Safety  System  Override  Switch  - SPST  toggle 

56  - Drive  Motor  Switch  - SPST  toggle 

57  - Water  Circulation  Pump  Switch  - SPST  toggle 

6 

58  - Oil  Circulation  Pump  Switch  - SPST  toggle 

59  - High  Pressure  Water  Pump  Switch  - SPST  toggle 

510  - Waviness  Drive  Switch  - SPST  toggle 

511  - Leak  Detector  Switch  - SPST  toggle 

REL1  - SPDT  Relay  in  CTL1 

REL2  - SPDT  Relay  in  CTL2 

REL3  - Safety  System  Relay  - SPDT  120  VAC 

REL4  - Leak  Detector  Relay  - SPDT  24  VAC 

PL1  - Pilot  Light 

PL2  - Pilot  Light 

Tl  - Leak  Detector  Transformer  - 24  VAC 

T2  - Cooling  water  servo  transformer  - 24  VAC 

FI  - Fuse 

CTL1  - Oil  Overtemperature  Control 

West  Model  10  0-300°F 

CTL2  - Seal  Face  Overtemperature  Control 

West  Model  10  0-500°F 

CTL3  - Cooling  Water  Temperature  Control  - Honeywell 

Dialatrol  R7352 
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Figure  4-13  (continued)  Machine  Electrical  Circuit 


Ml 

M2 

M3 

M4 

M5 

M6 

M7 

S2 

S4 

VI 

V2 

L1 

REL5 

REL6 

REL7 

TCI 

TC2 

TC3 


I 

- Drive  Motor,  Reeves  Vari-speed  Motodrive  Model  400 

- Water  Circulation  Pump  Motor,  Emmerson  JP540,  1 HP, 

3450  RPM 


- Oil  Pump  Drive  Motor,  Dayton  1725  RPM,  1/6  HP 

- Cooling  Fan  Drive  Motor,  1/40  HP,  1470  RPM 

- Waviness  Drive  Motor,  Dayton  Gear  Motor  6K301,  1/15  HP, 
30  RPM 


- Cooling  Water  Control  Servo  Motor,  Honeywell  M944B, 
Modutrol  Motor 

- Sump  Pump  Drive  Motor  AMXCO,  1/60  HP,  3000  RPM 

- Oil  Under  Pressure  Switch,  Honeywell  PA404B,  1023  2 

- Water  Pressure  Control  Pressure  Switch,  Honeywell  P634A 

- High  Pressure  Water  Pump  Electric  Solenoid  Air  Valve, 
Dayton  6X080,  1/4" 

- Leak  Detector  Water  Shutoff  Solenoid  Water  Valve,  Dayton, 
6X082,  1/2" 

- Leak  Detector.  Electrodes  insulated  by  dry  salt  soaked 
paper 

- Drive  Motor  Starter,  Dayton 

- Water  Pump  Motor  Starter,  Dayton 

- Oil  Pump  and  Fan  Motor  Starter,  Dayton 

- Thermocouple  in  Oil  Return  Line  - Type  J,  Omega 

- Thermocouple  Beneath  Seal  Face  - Type  J,  Omega 

- Thermocouple  in  Water  inside  Pressure  Vessel,  Type  J, 
Omega 


i 
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removed,  and  all  plumbing  and  hoses  hooked  to  test  ves- 
sel, balance  the  bridge  to  find  zero  torque.  Lock  all 
settings.  Using  the  calibration  stand,  check  range  cali- 
bration and  lock  all  settings.  The  torque  transducer  is 
to  be  left  on  for  the  entire  duration  of  the  test. 

2.  Lap  and  polish  the  carbon  and  secondary  ring  faces. 

Check  flatness  using  an  optical  flat.  Flatness  should 
be  within  one  light  band  across  the  width  of  the  face. 
Waviness  should  be  within  two  light  bands. 

3.  Make  total  roughness  and  CLA  traces  of  the  lapped  and 
polished  faces. 

4.  Mark  six  equal  spaces  around  the  carbon  ring  and  number 
one  through  six.  Using  an  appropriate  diameter  dowel 
pin,  determine  nose  height  to  0.0001"  or  better. 

5.  Install  seal  rings  making  certain  that  all  O-rings  are 
in  place  and  greased.  Wipe  face  of  seal  and  adjacent 
area  completely  free  of  grease  before  final  assembly. 

Make  certain  that  proper  shaft  spacer  is  used  to  com- 
pensate for  wear. 

6.  Install  thermocouples  making  certain  that  they  are 
clamped  in  position  so  as  not  to  vibrate  and  wear . 

7.  With  seal  dry  and  spring  load  only,  measure  and  record 
the  friction  torque  using  the  hand  crank  to  turn  the 
machine . 

8.  Install  bolt-tightening  bracket.  Install  pressure  ves- 
sel and  tighten  bolts  to  100  ft.  lb.  Remove  bolt-tight- 
ening  brackets. 


9.  Charge  accumulator  to  75%  of  operating  pressure  and  then 
close  accumulator  valve,  being  careful  not  to  overtighten 
valve . 

10.  Turn  on  main  central  switch.  Start  oil  pump.  This  will 
allow  other  equipment  to  be  operated.  Turn  on  all  water 
valves.  Water  valve  to  high-pressure  pump  should  only 
be  cracked.  All  other  valves  should  be  full  on.  Set 
oil  temperature  trip  to  160°F  and  seal  temperature  trip 
at  50°F  above  water  temperature.  Set  water  temperature. 

11.  Crack  two  bleed  valves,  one  below  and  one  above.  Start 
water  circulation  pimp.  Allow  water  to  flush  system 
completely  until  all  air  is  driven  out.  Turn  off  lower 
bleed  drain  valve  and  the  upper  bleed  valve . System  is 
now  ready  for  pressure.  Do  not  run  water  circulation 
pump  dry. 

12.  Turn  on  high-pressure  pump.  System  pressure  will  come  up. 

13.  Using  hand  crank,  measure  and  record  seal  torque. 

14.  Fill  leak  cavity  with  water  using  a small  diameter 
flexible  tube. 

15.  Apply  waviness  pressure  if  any. 

If.  Measure  static  leakage  if  any. 

17.  Start  chart  recorders,  and  check  pens  and  paper  supply. 

18.  Start  drive  motor,  bring  up  to  operating  speed  and 
check  speed. 

19.  Mark  time  zero  on  charts  and  check  leakage  on  graduated 
cylinder. 


20. 


Observe  that  system  temperature  stabilizes. 

21.  Record  leakage  as  necessary.  Record  steady  state  oil 
temperature  to  be  used  for  torque  zero  shift  corrections . 

22.  Run  machine  for  given  test  duration 

23.  After  test  is  completed,  reduce  drive  speed  to  low  and 
turn  off  all  systems  except  torque  transducer.  Drain 
pressure  vessel  and  system.  Mount  bolt-tightening 
bracket.  Remove  vessel  and  seal  assembly.  Note  wear 
condition  of  seal.  1 

24.  Remove  bolt-tightening  bracket  and  note  zero  torque 
position.  Compensate  for  temperature  if  any. 

25.  Using  previously  outlined  technique,  measure  nose  height 
of  carbon  seal  ring.  Measure  wear  on  secondary  ring  by 
making  several  traces  across  entire  width  of  face . 

26.  Make  total  roughness  and  CLA  measurements  in  radial  and 
tangential  directions  for  both  rings.  Make  traces  dia- 
metrically across  both  rings. 

27.  Correct  and  average  torque  data.  Compile  all  data  for 
test. 

EXPERIMENTS 

Several  different  series  of  experiments  have  been  planned 
using  the  test  apparatus. 

1.  Verify  the  waviness  model  by  experiment. 

Using  the  seal  test  rig,  experiments  will  be  conducted  to 
confirm  the  basic  predictions  made  by  the  model.  First,  the  effect 


of  waviness  amplitude  on  friction  and  leakage  will  be  determined. 
The  effect  of  waviness  on  wear  rate  will  be  established  by  com- 
paring the  wear  rate  for  a wavy  seal  to  the  wear  rate  for  a flat 
faced  seal.  Hydrodynamic  seal  performance  as  a function  of  pres- 
sure and  speed  will  also  be  established.  The  effect  of  surface 
roughness  on  hydrodynamic  lubrication  will  be  established  by  modi- 
fying the  natural  surface  roughness  of  the  mating  faces.  This  can 
be  achieved  by  introducing  abrasives  into  the  seal  system.  Surface 
characterization  will  be  performed  using  a surface  analyzer. 

2.  Demonstrate  the  traveling  wave  design  concept. 

In  order  to  take  advantage  of  waviness  to  reduce  friction 
and  wear,  it  is  necessary  to  move  the  wave  relative  to  the  seal 
face.  The  test  rig  seal  has  this  capability,  so  the  results  from 
(1.)  will  also  serve  to  demonstrate  the  advantages  and  practicality 
of  this  concept. 

3.  Experimentally  verify  a hydrostatic  seal  model. 

Reference  [23]  presents  a mixed  friction  hydrostatic  seal 

model  which  includes  radial  taper  and  wear.  This  model  has  been 
shown  to  theoretically  predict  important  seal  behavior  on  a quali- 
tative basis.  Thus,  it  is  thought  that  this  model  will  be  very 
useful  in  predicting  the  performance  of  non  wavy  face  seals.  In 
order  to  become  accepted  as  a useful  tool,  the  model  needs  to  be 
subjected  to  experimental  testing. 

The  seal  test  rig  will  be  modified  so  as  to  produce  a con- 
trollable radial  taper.  This  is  easily  accomplished,  given  the 
design  of  the  test  rig.  With  this  feature  it  will  be  possible  to 


be  presented  in  graphical  form  or  printed  form  for  further  inter- 


pretation. 

The  above  data  acquisition  system  will  also  be  used  to  ana- 
lyze the  surface  roughness  data  produced  by  the  surface  analyzer. 
Statistical  characteristics  of  the  surface  are  needed  in  conjunction 
with  studies  of  the  effect  of  roughness  and  questions  related  to 
roughness  modeling. 
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CONCLUSIONS 


Based  upon  the  work  presented  in  Chapter  2,  the  Christensen 
[20]  and  Tonder  [21]  model  for  isotropic  roughness  leads  to  conser- 
vative predictions  of  hydrodynamic  load  support.  Thus,  the  results 
given  in  Reference  [1] , the  previous  annual  report,  theoretically 
underestimate  the  percent  hydrodynamic  load  support,  and  can  be  con- 
sidered quite  useful  for  predicting  a lower  bound  on  hydrodynamic 
effects  in  the  presence  of  roughness. 

A truncated  roughness  model  was  studied.  By  making  certain 
assumptions,  this  model  avoids  certain  fundamental  questions  that 
have  been  raised,  while  still  providing  a theoretically  sound  model 
for  roughness  effects.  Results  based  upon  this  model  show  a higher 
percentage  of  hydrodynamic  load  support  than  the  previous  model. 
However,  experimental  data  is  needed  as  to  the  proper  truncation 
level  to  be  used  with  this  model.  This  data  will  be  supplied  by  a 
testing  program  currently  underway. 

In  Chapter  3 a design  procedure  for  a contacting  hydrodynamic 
face  seal  was  developed.  At  lower  numbers  of  waves  it  was  shown  that 
stability  problems  may  be  encountered,  while  at  higher  numbers  of 
waves,  these  are  limitations  imposed  due  to  the  physical  impossibility 
of  imposing  a large  enough  waviness.  While  this  procedure  is  not 
final  and  will  require  modifications  and  further  theoretical  and 
experimental  progress  is  made,  the  method  clearly  shows  the  details  of 
how  design  parameter  values  are  selected  to  produce  the  theoretical 
performance  obtainable  using  the  hydrodynamic  contacting  seal  concept. 
In  the  example  studied , a w ear  reduction  of  17:1  compared  to  a flat 
face  seal  was  shown  possible. 
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One  of  the  important  questions  remaining  is  whether  or  not 
sufficiently  small  operating  surface  roughness  can  be  obtained  to 
permit  the  type  of  hydrodynamic  operation  described.  This  question 
requires  an  experimental  answer. 

Finally  in  Chapter  4,  the  eiqjerimental  test  apparatus  is  des- 
cribed in  detail.  This  apparatus  will  be  used  to  study  the  hydrody- 
namic effects  above  as  well  as  certain  hydrostatic  effects.  The 
apparatus  operates  as  designed  and  is  currently  being  used  in  a test 
program  to  evaluate  the  effects  of  waviness  on  seal  performance. 


PUBLICATIONS 


Several  publications  have  resulted  from  this  research 
program  during  the  past  three  years.  Two  annual  reports,  ref- 
erences [1]  and  [5]  were  produced.  Two  papers,  [26],  [27]  based 
upon  the  first  one  and  one  half  year's  work  were  presented  and 
published.  Based  primarily  upon  work  performed  during  the  cur- 
rent contract  year,  two  Master's  theses  [28],  [29]  have  been  com- 
pleted. These  results  have  been  presented  as  parts  of  Chapters 
2 and  3.  A paper  is  being  written  based  upon  one  of  these  theses 
[28] . A paper  is  also  being  written  based  upon  the  results  in  the 
second  annual  report  [1] . 

During  this  writer's  Sabbatical  year,  two  papers  [23],  [24] 
were  written  and  submitted  on  certain  hydrostatic  effects  in  face 
seals.  Although  the  research  for  these  papers  was  not  performed 
as  a part  of  this  contract,  the  subject  matter  does  stem  from  the 


contract  research,  and  can  be  considered  as  a contribution  to  the 
total  effort  on  seal  research. 


During  September  1978  a presentation  on  seal  theory  was  made 
at  the  NASA  Lewis  Seals  Workshop. 
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